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ABSTRACT 


The  purpose  of  Section  I  Is  to  generalize  Theorem  3  of  Severo, 
Montzingo,  and  Schillo,  "Characterization  of  the  asymptotic  distribu¬ 
tions  of  a  transformed  normal  random  variable,"  Sankhya,  Series  A,  27, 
by  relaxing  the  assumption  of  normality  and  removing  the  requirement 
that  the  parameters  ^  and  of  the  distribution  of  the  random 

variable  Y  be  the  mean  and  variance,  respectively,  A  result  analogous 
to  the  above-mentioned  theorem  is  obtained  for  the  class  of  location- 
scale  parameter  distributions.  Examples  are  given  which  show  that  the 
conditions  given  are  sufficient,  but  not  necessary,  for  the  existence 
of  an  asymptotic  distribution  of  a  transformed  random  variable. 

Section  II  illustrates  ways  of  deciding  whether  or  not  a  given  uni¬ 
variate  random  variable  X  can  be  transformed  into  a  given  univariate 
random  variable  Y  ;  and  it  gives  procedures  for  defining  various  trans¬ 
formations  of  X  into  Y  in  the  event  that  one  such  transformation  is 
known  to  exist. 

Section  III,  which  consists  of  three  parts,  gives  an  illustration 
of  the  notion  of  robustness  of  a  test,  a  generalization  of  this  notion, 
and  a  tentative  definition  of  the  robustness  of  a  test  in  terms  of  a 
metric  on  the  space  of  power  functions  of  the  test.  A  short  investiga¬ 
tion  is  also  made  in  this  section  of  the  properties  of  the  Kolmogorov 
metric  on  the  space  of  location  families  of  distribution  functions,  and 
applications  are  made  '.o  the  normal,  Cauchy,  student' s-t,  gamma,  and 
exponential  distributions. 


ill 


Section  XV  presents  two  theorem*  that  provide  simple  interative 
solutions  of  special  systems  of  differentisl-difference  equations. 

Tha  first  system  consists  of  linesr  differentisl  equations  whose 
coefficient  matrices  are  triangular,  have  constant  elements,  and  have 
diagonal  elements  equal  to  each  other  at  most  in  pairs.  Tha  equations 
of  the  second  system  also  have  constant  triangular  coefficient  matrices, 
such  that  whenever  there  are  equal  diagonal  elements  then  sufficient  con¬ 
ditions  are  imposed  on  the  matrices  themselves  so  that  the  solutions 
involve  only  sums  of  exponential  terms. 

The  theorems  are  applied  to  the  simple  icochastic  epidemic  and  to  i 

the  general  stochastic  epidemic,  respectively,  in  each  of  which  the 
initial  distribution  of  the  number  of  uninfected  susceptlbles  and  the 
number  of  infective*  are  arbitrary  but  the  total  population  sire  is 
assuamd  bounded.  The  results  for  the  simple  stochastic  epidemic  provide 
solutions  not  obtainable  by  previously  known  results.  The  results  for 
the  general  stochastic  epidemic  are  simpler  and  more  direct  than  other  known 
methods,  which,  when  used  to  solve  the  problem  having  an  arbitrary  Initial 
distribution,  would  Involve  additional  steps  that  would  sum  proportionally- 
weighted  conditional  results. 
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I.  ON  THE  ASYMPTOTIC  DISTRIBUTION  OF  A  TRANSFORMED  RANDOM  VARIABLE 


In  [2],  there  la  stated  and  proved  a  theorem  regarding  the  asymptotic 
distribution  of  a  transformed  normal  random  variable,  and  the  following 
classes  of  Baire  functions  are  defined: 

For  each  natural  number  n  ,  and  each  real  number  s  ,  A  fn;a] 

is  that  class  of  Baire  functions  whose  elements,  X  ,  satisfy  the  con¬ 
ditions:  (1)  the  derivative,  X^  ,  of  X  is  continuous  at 

a  ,  (11)  X^n)(s)  is  a  nonzero  real  number,  and  (ill)  If  m  Is  a 

natural  number  less  than  n  ,  then  X'  (s)  Is  zero. 

For  each  natural  number  n  ,  and  each  real  number  s  ,  U[n;s] 

la  that  class  of  Baire  functions  whose  elements,  tu  ,  satisfy 

lim  tu(x)/xn  •  s/n! 
x-*0 

The  theorem  In  f 2 ]  to  which  we  refer  may  be  stated  In  a  slightly 
altered  form  as  follows: 

If  Y  Is  a  normal  random  variable  with  mean  u  and  standard 

y 

deviation  cr •  ;  if  Z  ■  (Y-|iy)/cry  ;  If  there  Is  a  natural  number  n 

such  that  h  £  A[n;ny]  ;  if  there  are  real  constants  R  and  Q  , 

with  Q  *  0  ,  such  that  r  £  u( n;Rh(n)  (yy)]  and  u0(cry)  -  b(uy)  +  r(cry)  ; 

If  £  0[n;Qh(n)(py)]  ;  If  X  ■  H(Y)  ,  and  If  W  -  fx-p0(ay)]/<r0(oy)  ; 

then  the  asymptotic  distribution  of  W  ,  as  cry  approaches  zero,  is 

that  of  the  random  variable  (Zn-R)/Q 

The  purpose  of  the  present  work  Is  to  relax  the  assumption  of 

normality  on  Y  and,  incidentally,  to  remove  the  requirement  that  the 

parameters  ^  and  7  be  the  mean  and  standard  deviation  of  Y  . 
y  y 
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For  convenience,  ve  introduce  sosae  terminology. 

Let  T  be  a  aubaet  of  Rj  X  r|  ,  the  real  upper  half  plane.  We 
call  a  family  (f(x;h,6) : (n,e)  £  T)  of  distribution  functions  a 
location-scale-parameter  family  with  respect  to  T  if,  and  only  if, 
there  is  a  distribution  function,  G  ,  on  the  real  line  such  that  for 
each  real  x  ,  ve  have  F(x;n,S)  »  G((x-p)/9)  ,  for  every  (n,0)  £  T  . 

For  ease  of  reference,  we  restate  Lemma  2  of  [2],  in  slightly 
altered  form. 

LEMMA.  Let  h  £  .AJnja]  *  R  be  a  real  number,  and 

Q  be  a  nonzero  real  number.  Then  a  necessary  and  sufficient 

condition  that  lim  [h(a+bz)  -  h(a)  -  r(b)]/q(b)  -  (zn-R)/Q  , 
b-0 

for  each  real  number  z  ,  is  that  r  £  Q[n;Rh^(a)]  and 
q  c  Qfn;Qh(n)(a)]  . 

We  can  now  prwe  a  generalized  version  of  Theorem  3  of  f  2 ] . 

THEOREM.  Let  Y  be  a  random  variable  with  distribution 
function  F(y;n,0)  ,  where  (F(y;u,9) : (n,&)  c  T)  is  a  location- 

scale  family  with  respect  to  T  .  Let  h  £  _A.fn;a]  ,  R 
be  any  real  number,  and  Q  be  any  n >nzero  real  number.  If 
W  -  f  h(Y)  -  h(y)  -  r(6)yq(0)  with  r  £  U[n;Rh(n) (u)]  and 
q  €  U[n;Qh^(y)]  ,  then  the  asymptotic  distribution  of  W  , 

as  9  approaches  zero,  is  equal  to  that  of..  (Zn-R)/Q  , 
where  Z  “  (Y-(i)/@  . 

Proof.  By  applying  Theorem  jS,  pag>  166,  of  [1],  the  characteristic 
function,  0^(0  »  “y  b«  written  in  the  form 

m  m 

PM(t)  -  J expt it  w(z)]dG(z)  -  j  exp{itfh(p+£z)  -  h(w)  -  r(0)]/q(6)}dG(x) 


where  G  ,  Che  distribution  function  of  Z  ,  doe*  not  depend  on  (u,6) 
By  applying  the  Lebesgue  dominated  convergence  theorem  and  the  Lemma, 
we  have,  for  each  real  t  , 

11m  P  (t)  -  /  exp[it(rn-R)/Q]dG(*)  -  3(t)  ,  »ay. 

6-*0  "  J 


Since  lim  3(t)  -  3(0)  ,  3  1*  continuous  at  zero  and,  hence,  3  is 

C-0 

a  characteristic  function.  Thus,  the  asymptotic  distribution  of  W  , 
as  8  approaches  zero.  Is  equal  to  the  distribution  of  (Zn-R)/Q  . 

We  note  in  passing  that  the  proof  does  not  depend  upon  the  existence 
of  moments  of  Y  .  For  example,  the  theorem  applies  to  the  situation  in 
which  Y  la  a  Cauchy  random  variable  with  location-scale  parameter  (p,6) 
belonging  to  T  *  Rj  X  r|  .  Thus,  for  any  transformation  h  of  the 
t>pe  described  in  the  thecrem,  the  asymptotic  distribution  of  the  random 
variable  W  »  fh(Y)  -  h(p)  -  r(<?)]/q (8)  is  equal  to  that  of  the  random 
/  variable  (Zn-R)/Q  ,  where  the  probability  density  function  of  Z  is 
given  by 


f(*) 


-*»  <  x  <  m  • 


We  also  note  chat  the  theorem  gives  conditions  under  which  the 
distribution  of  the  random  variable  W  approaches,  as  8  approaches 
zero,  the  distribution  of  the  random  variable  (Zn-R)/Q  which  does  not 
depend  upon  8  . 

We  give  below  some  examples  of  situations  in  which  the  distribution 
of  (Zn-R)/Q  may  depend  upen  8  ,  and  for  which  it  is  true  Chat  the 

distribution  functions  of  W  and  of  (Z°-R)/Q  approach  the  same  function 
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as  6  approaches  aero. 

gxanyle  1.  in  this  example,  v«  show  that  the  distribution  func¬ 
tions  of  V  and  of  (Zn-R)/Q  may  approach  the  same  distribution 
function  as  0  approaches  sero,  even  though  F(y;n»0)  is  not  a 
location-scale  family  with  respect  to  T  • 

Let  be  a  nonzero  real  number  end  let  0  be  a  real  number  in 
the  open  unit  Interval  (0,1)  •  Let  t'  ■  *Rd  let  Y(ii»0) 

denote  a  discrete  random  variable  with  distribution  function 


r(y;n.0)  • 


0  ,  -•  <  y  < 

02/2  ,  n-0  <  y  <  n 

l-02/2  ,  n  <  y  <  u+0 

1  ,  u+«  <  y  <  » 


If  ZCii»0)  denotes  the  random  variable  (Y(n,0)-n)/0  *  th<n  the 

tribution  function  of  Z(y,0)  is  given  by 


G(*;0)  - 


0  ,  —  <  *  <  -1 

e2/2  ,  -l  <  z  <  o 

i-e2/2  ,  o  <  z  <  l 

1  ,  1  <  z  <  • 


which  depends  on  (M,0)  ,  so  that  (F(y;M,0)  :  (H,0)  €  T’  }  is  not  a 

location-scale  family  with  respect  to  T' 

If  R  is  any  real  number  and  Q  is  any  positive  real  number,  then 

2 

we  denote  by  U(0)  the  random  variable  (Z  (y,0)-R)/Q  in  order  to 
stress  the  fact  that  the  distribution  function  of  depends  on 

6  .  It  is  easy  to  show  that  the  distribution  function  of  U(0)  is 


given  by 
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Thu*,  «a  6  approach**  s*ro,  th«  diatribucion  function*  of  the  random 
variable*  V(8)  and  U($)  both  approach  the  aame  distribution  function. 

Example  2.  In  this  example,  we  show  that  even  though  the  limit, 
a*  8  approaches  re rc\  of  the  distribution  function  of  the  random  variable 
U<0)  -  (Zn<M.S)-*)/Q  1*  not  a  distribution  function,  the  limit,  a*  8 

approaches  sero,  of  the  difference  betveen  the  distribution  function  of 
U(6)  and  that  of  V(6)  can  approach  sero. 

Let  T*  be  the  same  as  in  example  1,  and  let  Y0*,6)  denote  a 
discrete  random  variable  with  distribution  function 


r(y;u.s>  - 


0  ,  -•  <  y  <  U-0 

(1/2X1-S2)  ,  M-e‘l  <  y  <  U 

2  "  -I 

11/3(14$  )  ,  M  <  y  <  U4$ 

1  ,  U4 8  1  <  y  <  m 


The  distribution  function  of  the  random  variable  (Y(u,0)-u)/G  , 

which  we  again  denote  by  Z(u,@)  ,  is  given  by 


G(*;S) 


0  ,  —  <  s  <  -8~2 

(i/3(i-e2)  ,  *e*2  <  *  <  o 

(1/2X1+32)  ,  o  <  x  <  e~2 
l  ,  e2  <  z  <  - 


f 


2 

and  that  of  the  random  variable  U(0)  *  (Z  (u,6)-R)/Q  ,  for  positive 

real  R  and  positive  real  Q  ,  is  given  by 

t 


H(u;0) 


0  ,  —  <  u  <  -R/Q 

<  Q2  ,  -R/Q  <  u  <  (6'4-R)/Q 
1  ,  ( 8  4-R)/Q  5  u  <  - 


6 


By  taking  h  ,  r  and  q  to  ba  the  sawe  functions  aa  they  were 
In  exawpla  1  and  defining  W(0)  aa  It  was  defined  there,  we  find  that 
W(0)  haa  distribution  function 
f 


J(w;0) 


< 


0  ,  -•  <  w  <  -R(l+0)/Q (1+20) 

o2  ,  -R  (1+0) /Q  (1+20)  <  w<  [0_4-R (1+0) ]/Q (1+20)  . 

1  .  1 0'4-R (1+0)  1/Q  (1+20)  <  v  <  « 


In  order  to  simplify  the  confutations  and  sake  the  dlacuaalon  eaaler 
to  follow,  we  now  apeclallze  to  the  caae  R  •  1  and  Q  -  1  .We  note 
that  no  eaaentlal  generality  la  loat  by  ao  doing.  Under  thla  restric¬ 
tion  we  have 


H(u;0) 


0  ,  -•  <  u  <  -1 

02  ,  -1  <  u  <  0'4-l  , 

1  ,0  4  - 1  <  u  <  • 


and 


J(»;0)  - 


0  ,  —  <  w  <  -[1+0/ (1+0) l'1 

<  62  ,  -t  1+0/ (1+0) ]‘l  <  w  <  [0"4- (1+0)]/ (1+20)  . 

1  ,  [0“4- (1+0)]/ (1+20)  <  w  <  • 


It  la  clear  that,  for  any  fixed  real  number  x  ,  the  absolute 
difference  |j(x;0)  -  H(x;0)|  can  be  made  arbitrarily  swell  by  choosing 
0  sufficiently  close  to  zero;  l.e.. 


1 lm  | J (x ; 0)  -H(x;0)|  ■  0  ,  -•  <  x  <  •  . 

0-0 


It  la  also  clear  that  If  x  la  any  fixed  real  nuwber,  then  H(x;0) 
can  be  wade  arbitrarily  swell  by  choosing  0  sufficiently  cloae  to 
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1  1 
"Z  '« 

zero  (for  6  <  (x+1)  *  ,  H(xjS)  <  (x+1)  );  i.«., 

lim  H(x;6)  -  0  ,  -»  <  x  <  •  , 

6-0 

which  la  not  a  distribution  function. 
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II.  SUPPLEMENT  TO  CHAPTER  II  OP  ARL  65-75 


1 .  Preliminary  Remarks 

Let  B  denote  the  set  of  all  Balre  functions  which  map  the  real 
line  R  into  R  ;  and  let  0  denote  that  subset  of  B  which  contains 
the  function  P  if,  and  only  if,  F  is  a  distribution  function.  If  F 
is  a  function  in  D  ,  if  P  is  the  probability  measure  which  Induces 
F  ,  so  that,  at  each  point  x  of  R  ,  the  value  of  F  la 


» 


I 

j 


F(x)  -  z  :  z  <  x 

Vw 

if  h  is  a  function  in  B  ,  and  if  G  is  that  function  in  D  whose 
value  at  each  point  y  of  R  is 

G(y)  *  P  -f  x  :  h(x)  <  y 

i. 

then  the  symbol  (F,h,C)  stance  for  the  assertion  that  the  Balre  functions 
F  ,  h  and  G  are  related  to  one  another  In  the  manner  hypothesized. 

For  each  ordered  pair  (F,G)  of  functions  In  D  ,  the  symbol  (F,*,G) 
denotes  the  subset  -j"  h  :  (F,h,G)  j"  of  B  ;  and  (F,G)  is  said  to  be 
*  compatible  pair  or  an  incompatible  pair  according  as  It  is  not  or  is 
true  that  the  set  (P,*,G)  Is  the  empty  set  0  .  Chapter  II  of  the 

interim  technical  report  ARL  65-75  ,  entitled  "Some  Properties  Of  Dis¬ 

tribution  Functions  And  Transformations  That  Induce  One  Another",  pro¬ 
vides  necessary  and  sufficient  conditions  under  which  an  ordered  pair 
(F,G)  of  distribution  functions  in  D  is  a  compatible  pair.  In  this 
supplement  to  chapter  II  of  ARL  65-75  ,  we  shall  carefully  examine  some 

special  ordered  pairs  of  distribution  functions  in  D  ;  we  shall  determine 
whether  or  not  they  are  compatible  pairs;  if  one  of  them,  (F,G)  ,  is  a 

-9- 


-  «*!» 


compatible  pair,  then  ve  shall  determine  at  least  one  function  h  in 
the  non-empty  set  (F,*,G)  ;  and.  In  the  final  section  of  this  supplement, 

we  shall  consider  ways  of  finding  various  functions  in  (F,*,G)  when 
(F,G)  is  a  compatible  pair.  We  would  like  our  illustrations  to  require 
—  and,  thereby,  to  justify  --  some  of  the  elaborate  details  in  the 
theory  concerning  compatible  pairs  of  functions  in  D  which  is  developed 
in  chapter  II  of  ARL  65-75. 

In  order  to  stake  this  supplement  to  chapter  II  of  ARL  65-73  nota- 
t tonally  independent  of  that  technical  report,  we  shall  define  here  those 
of  its  symbols  and  terms  that  we  use.  Of  course,  many  statements  in  this 

V 

supplement  will  not  be  independent  of  that  reference,  because  we  shall 
not  prove  them  here. 

2.  The  Motion  of  Clearance 

Let  W  denote  the  set  of  all  infinite  column  matrices  whose  entries 

are  real  numbers.  Thus,  W  may  also  be  regarded  as  the  set  of  all  infinite 

sequences  of  real  numbers.  If  w  denotes  a  matrix  in  W  (or  a  sequence 

in  W  ),  then,  tor  each  positive  integer  n  ,  the  symbol  w^  denotes 

the  entty  in  the  n'th  row  of  w  (or  the  n'th  term  of  w  ). 

Three  operators  which  may  be  applied  to  symbols  denoting  elements 

of  W  are  cum:  ,  lim:  and  sum:  .  If  w  denotes  an  element  of  W  , 

then  cum:w  denotes  the  sequence  of  partial  sums  of  the  terms  of  the 

infinite  sequence  w  ;  that  is,  cum:w  denotes  that  element  v  of  W 

whose  first  term  v.  is  w,  ,  and  whose  n'th  term  v  ,  for  each  integer 
IX  n 

n  >  1  ,  is  wq  +  vn_j  .  If  w  denotes  a  convergent  sequence  in  W  , 

then  lim:w  denotes  the  limit  of  w  .  Finally,  if  v  denotes  an  element 
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of  W  which  is  such  that  cum:v  is  a  convergent  aequence,  then  sum:w 

denotes  the  real  number  lim:cua:v  . 

A  matrix  v  in  W  is  said  to  dominate  a  matrix  w  in  W  if,  and 

only  if,  for  each  positive  integer  n  ,  v  >  w  ;  and,  the  fact  that 

n  “  n 

v  dominates  w  may  be  asserted  symbolically  by  either  v»  w  or  w  «  v  . 

If  v  >>  w  and  v  +  w  ,  then  we  may  write  v  »  w  and  w  «  v  . 

Let  T  denote  the  set  of  all  matrices  with  infinitely  many  rows  and 

columns,  in  each  column  of  which  there  is  one,  and  only  one,  non-zero 

entry,  the  real  number  1  .  If  A  denotes  a  matrix  in  T  ,  then,  for 

each  ordered  pair  m,n  of  positive  integers,  the  symbol  A  denotes 

in  i  bl 

the  entry  in  the  m'th  row  and  n'th  column  of  A  .  Let  I  denote 

the  Identity  matrix  in  T  ;  that  is,  let  1  be  that  matrix  in  T  which  is 

such  that,  for  each  positive  Integer  n  ,  I  «  1 

— n,n 

For  each  matrix  w  in  W  ,  let  T(w)  denote  that  subset  of  T 

which  contains  the  matrix  A  in  T  if,  and  only  if,  the  matrix  product 

Aw  is  a  defined  column  matrix;  in  other  words,  let  T(w)  denote  that 

subset  of  T  which  contains  the  matrix  A  in  T  if,  and  only  if,  for 

each  positive  Integer  m  ,  the  sequence  v  in  W  ,  whose  n'th  term 

v  is  A  w  for  each  positive  integer  n  ,  is  such  that  sum:v  is 
n  id  ,  7i  n 

defined.  It  is  evident  that,  for  each  matrix  w  in  W  ,  T(w)  ^  0  , 

because  Iw  ■  w  . 

A  matrix  v  in  W  is  said  to  clear  a  matrix  w  in  W  if,  and 
only  if,  there  is  at  least  one  matrix  A  in  T(w)  for  which  it  is  true 
that  v  dominates  the  matrix  product  Aw  ;  and,  the  fact  that  v  clears 
w  may  be  asserted  symbolically  by  either  v  >»  w  or  w  «<  v 
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Our  special  interest  In  the  notion  of  clearance  --  that  la,  the 
idea  which  underliea  the  aaaertion  that  one  lnfinice  aequence  of  real 
mrabera  clears  another  such  aequence  —  will  be  reatricted  to  auch  a 
node at  range  of  ita  application  (e.g.,  to  auch  sequences  in  W  as  w  , 
where  w  is  a  monotone  non-increasing  sequence  of  non-negative  real  num¬ 
bers,  and  where  auai:w  is  a  definite  non-negative  real  number  that  does 
not  exceed  1  )  that  it  might  seem  that  some  simpler  notion  would  serve 

our  purpose  for  introducing  it.  It  happens,  however,  that  we  have  found 
no  replacement  for  it  that  is  as  easily  described  as  it  is,  and  that  is 
suitable  to  our  needs. 

3.  The  Theorem  in  Chapter  II  of  ARL  65-75 

Some  notation  from  chapter  II  of  ARL  65-75  is  the  following:  For 

each  real  number  x  ,  a  (x)  denotes  that  subset  of  R  which  contains 

the  point  z  of  R  if,  and  only  if,  z  <  x  .  For  each  function  F  in 

D  ,  the  symbol  R(F;»)  denotes  the  set  of  all  points  of  discontinuity 

of  F  ;  the  symbol  R'(F;°°)  denotes  the  difference  set  R  -  R(F;°*)  ; 

the  symbol  R  (F;«)  denotes  Chat  subset  of  R'(T;en)  which  contains  the 
c 

point  x  of  R'(F;<r)  if,  and  only  if,  for  each  point  z  of  R  which 

is  leas  than  x  ,  F(z)  <  F(x)  ;  the  symbol  R'(F;oc)  denotes  the  difference 

c 

set  R^F;00)  -  R^(?;er)  ;  the  symbol  Rw(F;»)  denotes  that  subset  of  R^CF;*-) 

which  contains  the  point  x  of  R  (F;*)  if,  and  only  if,  for  each  real 

c 

number  z  <  x  ,  and  with  P  denoting  the  probability  measure  that  Induces  F  , 
P  {  r  :  r  «  (0  (z)HR  (F;<r)  V  <  P  -j"  r  :  r  e  r»  (x)  r\  Rc(F;e n) 

^  J  V 

and  the  symbol  R'(F;«)  denotes  the  difference  set  R'(F;on)  -  R  (F ;») 

v  w 
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Finally,  for  the  function  F  in  D  and  the  real  number  x  ,  the  symbols 
R(F;x)  ,  R’(F;x)  ,  Rc<F;x)  ,  R^(F;x)  ,  Rw(F;x)  and  R^(F;x)  denote 

the  set  intersections  oj(x)n  R(F;»)  ,  a>(x)H  R’  (F;»)  ,  <u(x)n  R£  (F;«)  , 

<u(x)n  R^(F;<*)  ,  co(x)0  R^fF;®)  and  co(x)n  R^fF;®)  ,  respectively. 

Let  p  be  that  function  which  maps  D  into  W  in  such  a  way  that 
its  value  at  each  function  F  in  D  is  that  sequence  p(F)  =  w  in  W 
which  satisfies  the  following  six  conditions: 

(1)  if  R(F;«)  =  <t  ,  then,  for  each  positive  integer  n  ,  w  =*  n  ; 

n 

(2)  if  R(F;®)  contains  exactly  m  points,  the  greatest  of  which 

is  r  ,  then,  for  each  integer  n>m  ,  w  =r+n-m  ; 

n 

(3)  if  a  is  a  point  of  R(F;®)  ,  then  there  is  one,  and  only  one, 

positive  integer  n  such  that  w  =  a  ; 

n 

(A)  if,  for  some  positive  integer  n  ,  w  is  not  a  point  of  R(F;®)  , 

n 

then  R(F;®)  is  not  an  infinite  set; 

(5)  if  a  and  b  are  points  of  R(F;<r)  such  that  the  saltus  of  F 

at  a  is  greater  than  the  saltjs  of  F  at  b  ,  then  there  exist 

positive  integers  m  and  n  such  that  m  <  n  ,  w^  =  a  and 

w  =  b  ;  and 
n 

(6)  if  a  and  b  are  points  of  R(F;<r)  such  that  a  <  b  ,  and 

such  that  the  saltus  of  F  at  a  is  the  same  as  the  s*»ltus  of 

F  at  b  ,  then  there  exist  positive  integers  m  and  n  such 

that  m  <  n  ,  w  =  a  and  w  =  b 
m  n 

Let  W  be  that  subset  of  W  which  contains  the  sequence  w  of  W 
if,  ard  only  if,  cura:w  is  a  convergent,  monotone  non-decreasing  sequence, 
w  is  a  monotone  non-increasing  sequence,  and  sum;w  does  not  exceed  1 
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If  P  it  any  function  in  D  ,  let  v  be  the  sequence  p(F)  in 

W  ;  for  etch  positive  integer  m  ,  let  be  that  sequence  in  W 

whose  n'th  tern,  for  etch  positive  integer  n  ,  is  v^  m  F(w  )  - 

n  m 

P(wB  -  1/n)  ;  end  let  the  symbol  ssl:F  denote  that  sequence  u  in  W' 

whose  m'th  tern,  for  each  positive  integer  to  ,  is  u  =»  limtv^ 

n 

Thus,  sal:  is  an  operator  which  is  applicable  to  each  function  F  in 
D  ;  and,  loosely  speaking,  sal:F  is  the  monotone  non-increasing  se¬ 
quence  whose  non-zero  terns  are  the  saltuses  of  F 

By  making  use  of  the  notation  described  in  this  supplement  to 
chapter  II  of  ARL  65-75  ,  the  theorem  of  that  chapter  can  be  stated  con¬ 
cisely  as  follows: 

For  F  and  G  in  D  ,  (F,*,G)  ^  0  if,  and  only  if,  sal:F  <«  sal:G 

Thus,  this  theorem  gives  necessary  and  sufficient  conditions  for  an  ordered 
pair  (F,G)  of  functions  in  D  to  be  a  compatible  pair.  However,  since 
these  conditions  appeal  to  the  somewhat  formidable  notion  of  clearance  for 
their  meanings,  it  behooves  us  to  show  that  this  theorem  is  necessarily 
preferable  to  the  tautological  assertion  that  (F,G)  is  a  compatible  pair 
if,  and  only  if,  (F,G)  is  a  compatible  pair.  In  the  next  section,  we 
intend  to  provide  the  theorem  with  a  modest  justification. 

4.  Clearance  in  W* 

In  this  section,  we  shall  state  and  prove  three  theorems  concern¬ 
ing  the  clearance  relation  »>  between  sequences  in  W'  .  These  theorems 
will  be  considered  in  the  next  section  where  we  shall  be  dealing  with 
some  special  compatible  and  incompatible  pairs  of  discrete  distribution 
functions  in  D  .  However,  we  would  not  want  the  theorems  of  this 
section  to  be  mistaken  for  fundamental  theorems  in  a  carefully  constructed 
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theory  of  clearance.  Such  a  study  would  properly  begin  on  foundations 
provided  by  a  good  understanding  of  the  set  W  ,  and  of  the  dominance 
relation  »  between  its  elements,  and  by  a  separate  study  of  the  set  T 
In  unifying  the  results  of  these  separate  investigations  in  order  to 
develop  a  theory  of  clearance,  one  of  the  fundamental  theorems  would  assert 
that  clearance,  like  dominance,  is  a  transitive  relation.  The  following 
theorems  are  too  limited  in  scope  to  be  regarded  as  fundamental  theorems. 

THEOREM  1.  If  the  sequence  v  in  W  dominates  the  sequence  w  in  W'  , 
Chen  v  clears  w 

PROOF:  Since  v  >>  w  ,  ^  «  T(w)  and  I\i  m  v  ,  v  X>  Iw  ;  so 

V  »>  w  . 

THEOREM  2.  The  sequence  v  in  W'  clears  the  sequence  w  in  W  only 
if  cum:v  dominates  cum:w 

PROOF:  We  shall  prove  this  theorem  in  the  following  way:  we  shall 
let  v  and  w  be  any  sequences  in  W*  which  are  such  that  the  sequence 
y  a»  cum:v  in  W  dee*  not  dominate  the  sequence  x  =  cum:w  in  W  ;  we 
shall  let  A  be  a  matrix  in  T(w)  ;  and  we  shall  let  u  be  the  matrix 
product  Aw  in  W*  ;  then  we  shall  show  that  A  cannot  meet  the  necessary 
requirements  for  v  to  dominate  u  ,  so  that  v  cannot  clear  w 

Since  y  does  not  dominate  x  ,  there  exists  one,  and  only  one, 
positive  Integer  k  such  that  y^  <  x^  ,  and  such  that,  if  n  is  a 
positive  Integer  leso  than  k  ,  then  y^  j  x^  .  Since  v  and  w  are 
monotone  non-increasing  sequences  of  non-negative  real  numbers,  it  follows 
that 
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(*)  for  any  positive  integers  m  end  n  such  that  o  >  k  >  n  , 

v  <  v 

m  n 

For  each  positive  integer  n  ,  let  r(A;n)  denote  that  positive 

integer  m  which  is  such  that  A  «  1  ;  and,  for  each  positive  Integer 

0,11 

m  ,  let  A^  denote  that  sequence  in  W  whose  n'th  term  A^m^  is 
A^  n  for  each  positive  Integer  n  .  Now,  if,  for  some  positive  integer 
n  <  k  ,  the  integer  m  =  r(A;n)  is  not  less  than  k  ,  then,  by  (*)  , 


u  =  )  A^m^  w.  >  A  ^m^w  m  w  >  v 
m  L  h  h-nn  n  m 

h=l 


9 


so  that  v  cannot  dominate  u  =  Aw  ;  therefore,  if  v  is  to  clear  w  ,  A 
must  be  such  a  matrix  in  T  that,  for  each  positive  Integer  n  <  k  , 
r(A;n)  <  k  .  If  A  meets  this  requirement,  and  if  *  denotes  the  sequence 
cum:u  in  W  ,  then 


and.  Since  >  y^  ,  >  y^  .  Consequently,  there  is  a  positive 

integer  n  <  k  such  that  u_  >  v  ,  so  that  v  does  not  dominate  u 

—  n  n 

Thus,  for  no  matrix  A  in  T(w)  ,  does  v  dominate  Aw  ;  hence,  v 
does  not  clear  w 


THEOREM  3.  The  sequence  v  in  W*  clears  the  sequence  w  in  V  only 
if  sum:v  >  sum:w 
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PROOF:  Let  y  ■  cum:v  and  x  ■>  cum:w  .  If  aum:v  <  sum:w  ,  then, 
because  v  and  w  are  monotone  non-increasing  sequences  of  non-negative 
real  numbers,  there  is  a  positive  integer  k  such  that  >  sum:v  ; 

hence,  x^>  y^  ,  so  that  cum:v  does  not  dominate  cum:w  and,  by  theorem 
2,  v  does  not  clear  v 


5.  Examples  of  Compatible  and  Incompatible  Pairs 


Let  Iq  denote  the  negative  half  of  the  real  line  R  ;  let 
denote  the  non-negative  half  of  R  ;  and,  for  each  positive  Integer  n  , 

!l  I 

let  I  ,  I*  ,  J  and  J'  denote  the  respective  intervals  [  n  -  1  ,  n  -  1/2 
(  n  -  1/2  ,  n  )  ,  [  -2n  ,  -2n  +  99/100  )  and  [  -2n  +  99/100  ,  -2n  +  2  )  , 

which  are  all  open  on  the  right. 

Let  F  be  that  function  in  0  whose  value  at  each  x  e  R  is  that 


real  number  F(x)  which  is  defined  as  follows:  if  x  <  IQ  ,  then  F(x) 

“  101/2(101-x)  ;  if  n  ia  a  positive  integer  and  x  c  I  »  then  F(x)  * 

n 

2 

(101x+2n  +200n+10100)/2(n+100) (n+101)  ;  and,  if  n  is  a  positive  Integer 

and  x  e  I'  ,  then  F(x)  =  (4n2+602n+20099)/4(n+100) (n+101)  . 

n 

Let  H  be  that  function  in  D  whose  value  at  each  x  e  R  is  that 


real  number  H(x)  which  ia  defined  as  follows:  if  x  e  I- 


then 


H(x)  =  50/(100-x)  ;  if  n  is  a  positive  integer  and  x  «  I  ,  then 

n 

H(x)  =  (50x+n2+99n44950)/(n+99) (n+100)  ;  and,  if  n  is  a  positive  integer 

and  x  €  r  ,  then  H(x)  =  (n2+149n+4925) /(n+99) (n+100)  . 

Let  G  be  that  function  in  D  whose  value  at  each  x  e  R  is  that 


real  number  G(x)  which  is  defined  as  follows:  if  x  <  ,  then  G(x)  ■ 

(x+l)/(x+2)  ;  if  n  is  a  positive  integer  and  x  e  Jn  ,  then  G(x)  * 

(x+4n) /4n(n+l)  ;  and,  if  n  is  a  positive  integer  and  x  «  J'  ,  then 

n 

G(x)  =  (200n+99)/400n(n+l)  . 
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We  Intend  to  find  out  whether  or  not  (F,H)  end  (F,G)  are  com¬ 
patible  pairs.  Let  sal:P  -  w  ,  sal:K  -  u  ,  aal:G  ■  v  ,  cum:w  -  w  , 
cua:u  ■  u  and  cum:v  *  v  ;  and  let  a  ,  b  ,  a  and  b  be  thoae  sequences 

In  W  whose  n'th  terns,  for  each  positive  integer  n  ,  are  «  ■  u  -  w 

n  n  n 

b  ■  v  -  w  a  ■  u  -  w  and  b  ■  v  -  w  ,  respectively, 
nnn’nnn  n  n  n  '  r  7 

We  shall  consider  the  pair  (F,H)  first.  Each  of  the  sets  ..(F;*») 
and  R(H;<»)  is  the  set  of  all  positive  integers;  in  fact,  p(F)  and 
p(H)  are  the  sequence  of  positive  Integers.  Furthermore,  for  each  posi¬ 
tive  Integer  n  ,  w  -  101/4(n+100) (n+101)  ,  u  =  25 /(n+99) (n+100)  , 

n  n 

a  -  ( 101 -n)/4 (n+99) (n+100) (n+101)  ,  w  -  n/4 (n+101)  ,  u  »  n/4 (n+100) 

n  n  n 

and  a  «  n/4(n+100) (n+101)  .  Thus,  sum:w  =  lim:w  =1/4  and  sum:u 

n 

=1/4  .  In  theorems  2  and  3  of  section  4  ,  there  were  given  necessary 

conditions  for  u  to  clear  w  ;  since  u  dominates  w  and  sum:w  does 
not  exceed  sum:u  ,  these  necessary  conditions  are  met;  and,  therefore, 

It  cannot  be  concluded  that  (F,H)  is  an  incompatible  pair  by  applying 
thoae  theorems.  In  theorem  1  of  section  4  ,  there  was  given  a  sufficient 
condition  for  u  to  clear  w  ;  since  a  ^  ^  0  »  «o  that  u^q2  <  W102 
u  does  not  dominate  w  ,  this  svfflcient  condition  Is  not  met;  and,  there¬ 
fore,  It  cannot  be  concluded  that  (F,H)  is  a  compatible  pair  by  applying 
that  theorem.  In  this  example,  the  simple  tests  of  section  4  are  of  no  use 
to  us.  We  have  to  find  some  other  way  of  determining  whether  or  not  the 
set  (F,*,H)  is  empty. 

We  designed  this  problem  to  illustrate  some  of  the  difficulties  which 
can  attend  the  use  of  the  clearance  criterion  in  testing  the  compatibility 
of  ordered  pairs  of  distribution  functions  in  D  .  However,  since  we 


-18- 


wanted  a  decialon  to  come  out  of  a  modest  amount  of  study,  we  made  it 
possible  to  show  that  (F,H)  Is  an  Incompatible  pair  (i.e.,  that  u 
does  not  clear  w  )  in  a  way  that  is  only  a  little  more  complicated 
than  it  would  have  been  if  the  tests  derived  from  the  theorems  of  section 
4  could  have  been  employed  successfully.  Of  course,  in  this  way,  we  do 
no  more  than  hint  that  there  might  be  ordered  pairs  of  functions  in  D  to 
which  neither  the  adjective  "compatible"  nor  the  adjective  "incompatible" 
can  be  applied.  Until  an  upper  bound  on  the  possible  number  of  simple 
decisions  necessary  to  conclude  that  a  given  ordered  pair  of  functions 
in  D  is  or  is  not  a  compatible  pair  is  definitely  established,  the 
statement  that  the  pair  has  to  be  either  compatible  or  Incompatible  will 
appear  to  be  quite  flimsy  from  a  critical  viewpoint.  We  must  admit  that 
there  is  much  that  la  vulnerable  to  severe  criticism  in  our  simple  dichotomy 
of  the  Cartesian  product  set  D  x  D  of  all  ordered  pairs  of  functions 
in  0  into  compatible  and  incompatible  pairs;  but  no  example  here  is 
subject  to  that  kind  of  criticism. 

In  showing  that  (F,H)  is  an  incompatible  pair,  our  procedure  will 
be  as  follows:  we  shall  try  to  find  a  matrix  A  in  T(w)  which  is  such 
that,  if  x  is  the  matrix  product  Aw  in  W’  ,  then  u  dominates  *  ; 
therefore,  when  we  have  shown  that  no  such  matrix  A  can  meet  all  the 
requirements  which  we  shall  find  that  we  must  impose  on  it,  we  shall  have 
shown  that  u  cannot  clear  w  ,  and  that  (F,H)  is  an  incompatible  pair. 

For  each  positive  integer  n  ,  let  r(A;n)  denote  that  positive 

Integer  m  for  which  Am  n  ”  1  •  For  each  positive  Integer  n  ,  wn>  un+^  * 

because  w  -  u  ,  —  l/4(n+100) (n+101)  ;  therefore,  if  r(A;n)  were  to 

n  n-rl 

be  an  integer  m  greater  than  n  ,  then  the  m'th  term  of  the  sequence 


■  *  ‘■'-WwWWi 
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t  ■»  Aw  would  b« 


In  which  event  u  would  not  dominate  z  ;  consequently,  A  must  meet  the 
requirement  that,  for  each  poaltlve  integer  n  ,  r(A;n)  <  n. 

Since  r(A;l)  <1  ,  r(A;l)  ■  1  .  Suppose  that  It  has  been  shown 

that,  for  each  positive  integer  n  which  is  less  than  or  equal  to  a 
particular  poaltlve  Integer  k  ,  r(A;n)  must  be  n  if  it  is  to  be  at  all 
posaible  for  u  to  dominate  z  =  Aw  .  We  wish  to  see  whether  or  not,  under 
these  conditions,  r(A;k+l)  ,  which  cannot  exceed  k+1  ,  can  be  a  positive 
Integer  m  which  is  less  than  k+1  .  For  such  an  Integer  m  ,  the  m'th 

term  of  z  would  be 


-.-z 


A  w  >  A  w 
m,t  t  —  m,m  m 


m,k+l 


-  a  ) 
m 


If  w  .  la  greater  than  a  ,  then  z  would  be  greater  than  u  ,  and  u 
would  not  dominate  z  .  Hence  r(A;k+l)  cannot  be  an  integer  m  less  than 
k+1  if  exceeds  a^  .  Since,  for  each  positive  integer  m  <  201  , 

a  -  a  .  =  (201-m)/2(m+99)  (nrt-100)  (m+101)  (to+102)  is  positive,  and  since, 
for  each  integer  m  >  101  *  am  =  (101-ra)/4(m+99) (nH-100) (m+101)  is  negative, 

the  greatest  terra  in  the  sequence  a  of  W  is  a^  =  1/41208  .  Therefore, 
since  w  must  not  exceed  a  <  a  if  it  is  to  remain  possible  for  m 

to  be  less  than  k+1  ,  it  follows  that  m  cannot  be  less  than  k+1  if 


^k+l  exceet*8  ai  •  We  find  that  the  least  positive  integer  k  for  which 
it  is  true  that  does  not  exceed  a^  is  the  least  positive  integer  k 
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such  Chat  101/4(k+101) (k+102)  <  1/41208  ;  thus,  we  conclude  that  k  must 


exceed  (V4162009  -  203)/2  ,  which  means  that  k  >  919  .  For  thia  reaaon,  I 

“  1 

we  are  forced  to  :onclude  that,  if  k  were  to  be  leaa  than  919  ,  then  • 

r (A ; k+1 )  would  have  to  be  k+1  ;  consequently ,  for  each  positive  integer 
n  <  920  ,  A  must  meet  the  requirement  that  r(A;n)  ■  n  .  However,  if  A 

meets  thia  requirement,  then,  aince  *2.02  <  ®  102'nd  term  of  z  *  Aw 

ia 

00 

*102  “  X  A102,t  wt  -  A102,102  W102  “  W102  “  U102  "  *102  >  u102  , 

t=l 

ao  that  u  does  not  dominate  z  .  Thus,  no  matrix  A  in  T(w)  ia  such 
that  u  dominates  Aw  ;  consequently,  u  does  not  clear  w  ;  and  (F,H) 
ia  an  incompatible  pair. 

Now,  let  ua  consider  the  ordered  pair  (F,G)  .  It  will  not  be  hard 

for  us  to  show  that  (F,G)  is  a  compatible  pair,  so  that  the  subset  (F,*,G) 

of  B  is  not  9  .  Our  work  on  this  compatible  pair  will  be  primarily  that 

of  defining  a  particular  Balre  function  h  in  (F,+,G)  .  But,  first,  we 

shall  take  the  trouble  to  show  that  sal:G  *  v  clears  sal:F  **=  w 

We  find  that  R(G;®>  is  the  set  of  all  non-positive  even  Integers, 

that  p(G)  is  the  strictly  monotone  decresslng  sequence  of  non-positive  even 

Integers,  and  that  the  n'th  terms  of  the  sequences  v  ,  v  ,  b  and  b  , 

for  each  positive  integer  n  ,  are  v  *  101/400n (n+1)  ,  v  **  101n/400(n+l)  , 

n  n 

b  *v  -w  =  101<10100-H01n-99n2)/400n(n+l)(n+100)(n+101)  and 
n  n  n 

b  “  v  -  w  =  n(n+1010iy400(n+l) (n+101)  ,  respectively,  so  that  sum:v  «* 

n  n  n 

llm:v  *  101/400  .  Evidently,  the  tests  derived  from  the  three  theorems  of 

section  4  fail  to  disclose  whether  or  not  <F,G)  is  a  compatible  pair. 
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V«  shall  show  that  v  clears  w  by  defining  a  matrix  A  In 

T(w)  which  la  such  that  v  dominates  Aw  ,  This  satisfactory  matrix  A 

will  then  be  of  use  to  us  In  our  work  defining  a  function  h  In  (F,*,G)  . 

For  each  ordered  pair  of  positive  Integers  (m,n)  ,  let  A  be  either  1 

m,n 

or  0  according  as  It  Is  or  la  not  true  that  lOl(m-l)  <  n  <  101m  .  If  s 

Is  the  matrix  product  Aw  In  W*  ,  then,  for  each  positive  Integer  m  , 
the  m’th  term  of  z  la 


consequently,  v  dominates  z  ,  v  clears  w  ,  and  (F,G)  Is  a 
compatible  pair. 

The  matrix  A  establishes  the  following  one-to-many  correspondence 

between  the  points  of  R(G;<*>)  and  R(F;»)  ;  to  each  positive  integer  m  , 

and,  hence,  to  the  m'th  term  y  »  2  -  2m  of  p(G)  ,  there  correspond 

m 

101  points  of  R(F;«>)  ,  the  n*th  of  which  is  x  =  101(m-l)  +  n  of 

m,n 

R(F;<*>)  for  each  positive  Integer  n  <  101  ;  and,  to  each  positive  Integer 
n  ,  and,  hence,  to  the  n'th  term  n  of  p(F)  ,  there  corresponds  one, 
and  only  one,  point  of  R(G;»)  ,  the  even  integer  2y  which  is  such  that 
(l-n)/101  <  y  <  (101-n)/101  .  We  are  now  prepared  to  begin  defining  a 
function  h  In  (F,*,G)  by  following  the  procedure  used  in  chapter  II  of 
ARL  65-75  . 


For  any  subsets  M  and  N  of  the  real  line  R  ,  let  the  symbol 
B[M,Nj  denote  that  subset  of  B  which  contains  the  Bdire  function  q  In 
B  if,  and  only  If,  t'ae  value  of  q  at  each  point  x  of  M  Is  a  point 
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q(x)  of  H  .  Let  y  denote  that  apeciel  function  In  B  whose  value  at 
each  point  x  of  R  la  that  integer  y(x)  which  la  such  that 
x  -  1  <  T  (x)  <  x  .  Since  the  function  h  in  (F,*,G)  which  we  want  to 
define  la  a  function  in  B[R(F;«°),R(G;*»)]  ,  we  make  uae  of  the  many -to- one 
correspondence  between  the  points  of  R(F;»)  and  R(G;«)  which  waa  defined 

in  the  last  paragraph;  and  we  define  the  value  of  h  at  each  point  x  of 

R(F;«)  to  be 

(1)  !>(«)  . 

Of  course,  if  we  were  to  set  up  a  different  many-to-one  correspondence 

between  the  points  of  R(F;»)  and  the  points  of  a  subset  of  R(G;<») 

(l.e. ,  if  we  were  to  find  a  matrix  C  A  in  T(w)  such  that  v  >>  Cw)  , 

then  we  would  have  to  define  h  differently  over  R(F;«)  . 

The  subsets  R'(F;«)  and  R'(F;«)  of  R  are  the  same  set;  this 
c  w 

set  contains  the  point  x  of  R  if,  and  only  if,  there  is  a  positive 

Integer  n  such  that  I'  contains  x  ,  and  such  that  n  is  not  x  +  1/2  . 

n 

For  each  point  x  of  R'(F;<»>)  =  R*  (F;«>)  ,  there  is  a  real  point  y  >  x 

c  w 

such  that  F(y)  *  F(x)  .  We  are  free  to  define  h  over  R^'F;*)  in  any 
way  that  we  find  convenient;  therefore,  we  define  the  value  of  h  at  each 
point  x  of  R^(F;»)  to  be 

(2)  h<„)  .  -2  r  (  )  . 

It  remains  for  us  to  define  h  over  Rw(F;»)  ■  -  R^(F;»)  . 

Let  p  be  that  sequence  in  W'  whose  n'th  term,  for  each  positive  integer 


23 


d  ,  is  p  «v  -  ■  ■  l/400n(n+l)  ;  then  the  specisl  number  t  ■  sum:p  is 
n  n  n 

1/400  .  Let  P  denote  the  probsbility  measure  which  induces  the  distri¬ 
bution  function  P  ;  and  let  T  and  S  be  Borel  sets  in  R  which  ar'  such 
that  THS-0  ,  TUS«Rw(F;«)  and  PfT}  -  t  .  Since  R  -  R(F;®)U  R'  (F;») 
and  R'CFjco)  ■  R^(F;»)URw(F;»)  -  R^(F;«)UTUS  .since  THS=0  and 
*(F;»)r>R^(F;«)  -  R(F;®)f>T  -  R(Fi«)AS  -  ^(F;»)OT  -  R^(F;®)r»  S  -  0  , 

and  since  PfRl  *»  1  ,  PfR(F;®)l  «  sum:w  ■  1/4  ,  PfR^(F;«)]  ®  0  ,  and 

P[T1  -  t  -  1/400  ,  it  follows  that  R  *  R(F;«)U^(F;»)UTUS  and  P[Rl  * 
PfR(F;®)]  +  P[R'(F;®)]  +  PfT]  +  PfS]  ,  so  that  the  special  number  s  =  Pfli] 

i  i 

is  299/400  .  There  are  non-denumerably  infinitely  many  ways  of  choosing 

the  disjoint  Borel  sets  T  and  S  so  that  TUS  ®  RW(F;®)  ,  PfTj  =  t  and 
PfSl  *■  s  ;  and  no  two  such  choices  would  yield  the  same  definition  of  h 
over  Rw(F;®)  .  We  shall  define  S  (and,  consequently,  T  )  in  a  way 

that  we  feel  is  most  convenient. 


Let  SQ  and  denote  the  set  of  all  non-positive  real  numbers.;  and, 

for  each  positive  Integer  m  ,  let  S  denote  the  open  sub- interval  (  m  -  1  , 

m 

m  -  101/200  )  in  I  ,  and  let 
m 


S’  = 
m 


U 


1=0 


Since,  for  any  Integers  m  and  n  such  that  m  >  n  >  0  ,  =  0  , 

SmC  Rw(p;®)  ar.d  SnC  Rv(F;n)  Q  Rw(F;m)  ,  it  follows  that,  for  these  integers, 

PfS  US  1  =  PfS  1  +  Pf S  1  and  S'  C  S '  C  R(F;m)  .  Therefore,  since  q,,  = 
cd  n  in  n  n  m  v 

Pf SqI  ■  F(0)  «  1/2  and,  for  each  positive  integer  m  ,  q^  =  PfS^I  =  F(m  - 
101/200)  -  F(m  -  1)  =  9999/400(m+100) (m+101)  ,  it  follows  that  q^  =  PfS^l 
=  Pf SqI  =  1/2  and,  for  ecch  positive  integer  m  , 
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.  ,  Til  1  V1  .f.  i  1  .  99m  299  9999 

<C  «  Pfs^]  -  P  [  St  J  ptsil"  2  +  400(m+101)  “  400  '  400(nrt-101)  * 

1-0  1-0 

The  fact*  th*t,  if  q*  Is  that  sequence  In  W'  whose  a'th  term  Is  q^ 

for  each  positive  Integer  m  ,  then  lim:q*  -  29V/400  ,  and  that  S^C  Rw(p»") 

for  each  non-negative  Integer  m  ,  combine  to  enable  us  to  define  the  set  S 

to  be  the  following  union  of  Borel  sets  in  Rw(F;®)  : 

00 

U  si  • 

i«0 

For  each  positive  integer  m  ,  let  T”  denote  the  closed 

OR 

interval  [m  -  101/200  ,  m  -  1/2]  in  R  ;  let 


iui  m 

T  -  IJt!'  ;  and  let  T'  =  U  T.  . 

m  101  (m-l)-H.  m  ^  i 

Since,  for  any  Integers  m  and  n  such  that  m  >  n  >  1  ,  T^nT”  **  9  > 

T  "f*  (F;m)  and  I"(R  (F;n)  C  R  <F;m)  ,  it  follows  that,  for  these 

n  v  n  w  w 

integers,  P^’  'U  T"  1  -  P[T’ '  1  +  PfT' '  1  ,  T  C  R(F;101m)  and  T  C  R  (F;101n)  . 

m  n  d  n  m  w  41  w 

Also,  for  any  integers  m  and  n  such  that  m  >  n  >  1  ,  T^OT^  *  0  »  80 

that  P[T  UT  1  *  P[T  1  +  PfT  ]  .  Therefore,  since,  for  each  positive  integer 
m  n  m  n 

m  ,  p”  -  P[T"1  -  F(m  -  1/2)  -  F(m  -  101/200)  =  101/400(m+100)  (m+101)  ,  it 
in  ® 

follows  that 

r  i°i  l  ,  i 

Pm  “  PtTm1  “  P  L  H  T101(m-1)-H  J  =  ?  C  T101  (m-D+i1  “  400m(m+l) 

i-1 

for  each  positive  integer  m  .  Furthermore,  for  each  positive  integer  m  , 

m  m 

-f[U  I,]  -  £  *1,1  "ife  -4ooW  ‘ 

1=1  i-1 

The  fact  that  S  is  the  set  of  all  non-positive  real  numbers,  and  that,  for 
0 


each  positive  integer  a  ,  S  (IT"  ■  0  ,  S  III"  «  («  •  1  ,  1/2]  , 

an  an 

which  is  (I  Ul'  ]nR*  <F;a-l/2)  ,  and  (  m  -  1/2  ,  m)  -  I'AR^F;®)  . 

an  a  c 

forces  us  to  define  T  to  be  the  following  unions  of  Borel  seta  in  Rw(F;«) 

0  Ti'  -  U  T  . 

i-1  1  o-l  ■ 

The  most  convenient  definition  of  h  over  T  is  the  following: 

if  x  is  s  point  of  T  ,  and  if  m  is  that  positive  Integer  such  that 

contains  x  ,  then  the  value  of  h  at  x  is  h(x)  =  y  =  2(l-m)  .  How- 

a 

aver,  other  definitions  are  possible;  for  example,  by  observing  that 

11  139 

£  p  fT*  M  »Y  p[t*  ' ]  , 

i=5  MlO 

and  by  letting 

11  139 

M  «  T"  ,  T"  ,  T*  «  (T.UN)  •  M  , 

1*6  1  1=110  1  1 


T2  "  (T2UM)  “  N  *nd> 

for  each 

Integer 

m 

>  2  , 

T*  =■  T  ,  we  could  define 

m  m 

h  over  T  as  follows: 

if  X 

is  a  point 

of  T 

,  and  if  m  is  that  posi- 

tive  integer  such  that 

T*  contains  x 

■ 

» 

then 

h(x)  =  y  .  Nevertheless, 

m 

our  choice  of  definition  of  h 

over  T 

has  the 

advantage  that  its  value 

at  each  point  x  of  T 

lc 

(3)  h(x)  =-2  r 

( 4x-l  N 
\  A  04  / 

)  • 

Since  the  formulas  (1)  ,  (2)  and  (3)  are  all  the  same,  the  function  h 
has,  thus  far,  been  defined  so  that  it  coincides  with  a  convenient  monotone 
non-increasing  step-funct  .on  over  the  Borel  set  R-  S  . 

It  remains  for  us  to  define  h  over  S  .  Such  a  definition  will  be 
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given  at  the  end  of  a  process  which  follows  that  In  the  proof  of  the  theorem 
In  Chapter  II  of  ARL  65-75  . 

First,  we  define  the  function  K  to  be  that  function  which  maps  R 

Into  the  set  of  all  Borel  subsets  of  the  set  TUR(F;»)  in  such  a  way 

that,  at  each  point  y  of  R  ,  the  value  of  K  is  the  Borel  set 

K(y)  =  (x  :  x  e  TUR(F;»)  and  h(x)  <  y)  . 

In  order  to  put  this  definition  of  K  into  a  more  tractable  form,  ws  define 

the  following  sets:  for  each  positive  integer  m  ,  let  X  he  that  subset 

01 

of  R(F;<»)  which  contains  exactly  101  points,  the  n'th  of  which  Is 

x  =  101(m-l)+n  ;  let 

m,n 

m  m  oo 

x ’  =  X  ;  let  L  «  T  U  Xm  ;  let  L •  -  L  ;  and  let  L  =  \J  L 

"  1=1  1  m  m  m  »  i_l  1  1 

Since  L'  »  T'lJX1  for  each  positive  Integer  m  ,  L  »  TU8(P;«)  .  Thus, 

si  m  bb 

If  y  is  a  point  of  JQ  ,  then  K(y)  =  L  ;  and,  If,  for  some  positive 

Integer  ■  ,  y  Is  a  point  of  J  UJ'  ™  [-2m  ,  2  -  2m)  ,  then  K(y)  «  L  — 

ib  m 

L’  .  Since  PfLl  =  101/400  and,  for  each  positive  Integer  m  ,  PfX  1  = 

tn  ib 

z  «  l/4m(m+l)  , 

LU 

■ 

prxM  -  £  pfxti  -  .»0  -  prT;i  +  pr*;i  -  ■ 

1=1 

It  follows  that.  If  y  e  JQ  ,  then  P[K(y)l  =  101/400  ;  and,  if,  for  some 

positive  Integer  m  ,  y  €  J  UJ'  •  then  P[R(y)1  =  101/400(m+l)  . 

ib  vi 

Next,  we  define  the  function  g  to  be  that  'unction  In  B  whose 
value  at  each  point  y  of  R  is  g(y)  =  G(y)  -  PfK(y)1  .  Thus,  if  y  e 
JQ  ,  then  g(y)  •  (299y+198)/400(y+2)  ;  and,  If,  for  some  positive  Integer 

m  ,  y  c  J  U  J1  ,  then  g(y)  Is  either  (100y+299m)/400m(»fl)  or  99/400u« 

m  m 
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according  as  y  la  or  la  not  leaa  than  99/100  -  2m  . 

Now,  we  Intend  to  define  a  function  V  which  mapa  the  closed  Interval 

[0  ,  al  ,  where  a  »  P[Sl  «  299/400  ,  into  the  set  of  all  Borel  subsets 
of  S  In  such  a  way  that  each  of  the  following  four  conditions  is  satis- 
fled; 

(1)  for  each  point  it  of  fO  ,  s]  ,  the  image  V(u)  of  u 
under  the  mapping  V  is  a  Borel  subset  of  S  such  that 
P[V(u)l  -  u  ; 

(ii)  for  any  points  u  and  v  of  [0  ,  si  such  that  u  <  v  , 

V(u)  is  a  proper  subset  of  V(v)  ; 

(iii)  if  u  *  0  ,  then  V(u)  =  0  ;  and 

(iv)  if  u  *=  s  ,  then  V(u)  =  S  . 

The  particular  function  V  which  we  shall  define  is  taken  from  a  non- 
denumerably  infinite  collection  of  functions,  each  of  which  r  tisfies  these 
four  conditions,  but  no  two  of  which  yield  the  same  definition  of  h  over 
S  . 

Let  X  denote  that  subset  of  R  -  S  which  contains  the  point  x 
of  R  -  S  if,  and  only  if,  there  exists  a  positive  integer  m  such  that 
x  *m  -  101/200  ;  and  let  y  be  that  function  which  maps  SUX  into  the 
class  of  all  Borel  subsets  of  S  in  such  a  way  that  the  image  of  each  point 
x  «  SUX  under  the  mapping  y  is  the  set  u(x)  =  SAR'(F;x)  if  x  €  S  , 
and  is  the  set  n(x)  =  if  x  is  the  point  m  -  101/200  of  X  .  Thus, 
for  each  point  x  of  SUX  ,  either  x  <  0  and  P[u(x)l  =  F(x)  ,  so  that 
?fu(x)l  -  101/2 (101-x)  ,  or  there  exists  a  positive  integer  m  such  that 

m-1  <  x  <  m  and  P[u(x)1  “  ^  ~  ^  +  2wffl(x-m+l)  =  +  2wn(x-nri-l)  ,  so 
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If  X 


that  Pru(x)l  -  (20200x+299m2+29900m+2030201)/400(n+100)(m+101)  . 

and  z  are  points  of  SUX  such  that  x  <  z  ,  then 

0  -  Pf 01  <  Pfu(x)1  <  Pfp(z)1  <  P[Sl  -  s  , 

so  that  u(x)  Is  a  proper  subset  of  p(z)  .  Consequently,  the  function  u 
is  a  one-to-one  napping  of  Dorn  p  ,  its  domain  of  definition  SUX  ,  onto 
Ran  p  ,  its  range  of  values,  which  is  a  proper  subset  of  the  set  of  all 
proper  Borel  subsets  of  S  ;  in  other  words,  p  defines  a  one-to-one 
correspondence  between  the  elements  (points)  of  Dom  p  and  the  elements 
(sets)  of  Ran  p  .  Furthermore,  there  is  a  one-to-one  correspondence 
between  the  points  of  the  open  interval  (0  ,  s)  and  the  points  of  the  set 
SUX  =  Dom  p  such  that  u  e  (0,  a)  and  x  e  SUX  correspond  to  one 

another  if,  and  only  if,  Pfp(x)1  =  u  .  Therefore,  there  is  a  one-to-one 

correspondence  between  the  points  of  (0,  s)  and  the  Borel  sets  in  Ran  p 
such  that  the  point  u  c  (0  ,  s)  and  the  set  Z  c  Ran  u  correspond  to  one 
another  if,  and  only  if,  PfZl  ■  u 

We  define  the  function  V  ,  whose  domain  of  definition  Dom  V  is 

the  closed  interval  [0  ,  a]  ,  in  such  a  way  that  its  range  of  values  Ran  V 

is  Ran  p  augmented  by  the  values  V(0)  =  0  and  V(s)  =  S  of  V  at  0 
and  s  ,  respectively,  and  its  value  at  each  point  u  of  (0  ,  s)  *  Ran  p 
is  that  sat  V(u)  =p(x)  in  Ran  p  ,  for  some  unique  point  x  of  Dom  u  , 

which  is  such  that  P[p(x)1  =  u  .  Thus,  if  u  «=  0  ,  then  V(u)  =»  0  ;  if 

0  <  u  <  j  ,  then  V(u)  =  u  f  ^  for  -0Be  positive  integer  m  , 

.  ^  .  .  . .  {  400(m+100)(m+101)u-299m2-29900m- 2030201  N 

Ci<u^<*then  v<u)=u 1 ^ — 20200 - ;  5 

and,  if  u  =  s  =  299/400  ,  then  V(u)  =  S  . 


Let  C  be  tbet  function  which  meps  R  onto  Ran  p  in  such  a 

way  that  ita  value  at  each  point  y  of  R  ia  C(y)  »  V(g(y))  .  Thua,  if, 

99 

for  aoiae  positive  integer  a  ,  -2a  <  y  <  -Jqq  -  2m  ,  then 


-m( 


100y+299a 


)  ;  if. 


for  some  positive  integer  m  , 


-2a<y<2-2a,  then  C(y)  =  y  ^ 


101(99-200nQ 

99 


)  ;  if  0  <  y  , 


then  C(y)  *  4  ^ 


299y+198 


^  ;  and,  if,  for  some  positive  integer  m  , 


400m+20002  .  „  400m+20402  ..  .  f 

- - <  y  < - — -  ,  then  C(y)  =  M  ^ 


299m49799  2(m  +201m+10100) 

200  "  101 (y+2) 


Let  C*  be  that  function  which  maps  R  into  the  set  of  all 

aubsets  of  R  in  such  a  way  that  its  value  at  each  point  y  of  R  is  that 

subset  C'(y)  of  C(y)  which  contains  the  point  x  of  C(y)  if,  and  only 

if,  for  each  real  number  z  <  y  ,  x  is  not  s  point  of  C(z)  .  Thus,  if,  for 

some  positive  integer  a  ,  -2a  <  y  <  2  -  2m  ,  then  C'(y)  is  either  the 

2 

empty  set  0  or  the  set  whose  only  point  is  101(100y-200m  +99m)/(100y+299m) 
according  as  y  does  or  dots  not  exceed  99/100  -  2m  ;  if  0  <  y  <  202/99  , 
then  C'(y)  is  the  set  whose  only  point  is  101(99y-202)/(299y+198)  ;  and, 
if,  for  some  positive  Integer  a  ,  (400m+20002)/9999  <  y  <  (400m+20402)/9999  , 

.  .  ...  299td+9799  2(a2+20lm+10100) 

then  C*  (y)  is  the  set  whose  only  point  is  - joo -  "  — \0l(y+2) -  ’ 

Finally,  we  define  the  function  h  over  the  set  S  in  such  a 
way  that  ita  value  at  each  point  x  of  S  is  that  point  h(x)  =  y  of  R 
which  is  such  that  C'(y)  «=  x  .  Thus,  if,  for  some  positive  integer  m  , 


101 (200m-frl01)  _  101 (200m- 99)  a(299x+20200m-9999)  . 

99  -  ^  99  »  tnen  lOO(lOl-x) 


30- 


“  <  X  <  0  ,  then  h(x)  **  ~999^299'x^'  ’  and*  if*  for  aome  P°*ltlve 


integer  m  ,  m 


1  <  x  <  m  - 


101 

200 


.  ,  .  400a01x-hn2+50m+5151)+2m+202 

th€n  H(X> -  101(299m+9799-200x) - 


There  can  be  little,  doubt  that  the  method  used  in  this  construction 
of  a  function  h  in  (f,*,G)  leaves  much  to  be  desired.  After  having  shown 
that  (F,G)  ia  a  compatible  pair,  we  might  have  found  an  easier  way  of 
defining  a  function  in  (F,*,G)  .  However,  it  does  not  seem  likely  that 

there  is  a  procedure  which  is  necessarily  always  most  convenient  for  dealing 
with  problems  of  this  kind.  In  fact,  this  particular  illustration  in¬ 
volves  very  few  of  the  difficulties  which  are  anticipated  in  the  theory  of 
Chapter  II  of  ARL  65-75  . 

Of  course,  as  was  mentioned  earlier  in  this  section,  it  is  possible 
that  a  decision  on  the  matter  of  whether  or  not  a  given  ordered  pair 
(F,G)  of  functions  in  D  is  a  compatible  pair  might  be  precluded  by 
some  inherent  difficulties  in  the  definitions  of  F  and  G  .  This  possi¬ 
bility  arises  when  one  takes  a  mathematical  intuitlonist' a  critical  point 
of  view  and  observes  the  shaky  foundations  under  the  measure-theoretic 
notions  which  were  entertained  in  the  proof  of  the  theorem  of  Chapter  II 
of  ARL  65-75  . 


Now,  once  it  has  been  shown  by  a  satisfactory  procedure  that  a  given 
ordered  pair  (F,G)  in  0  x  D  is  a  compatible  pair,  and  that  a  definite 
function  h  is  in  the  non-empty  set  (F,*,G)  ,  the  matter  of  defining 
other  functions  than  h  in  (F,*,G)  can  be  considered.  Th«  next  section 
ia  devoted  to  this  consideration. 
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6.  Some  Members  of  the  Class  C ^ 

The  class  is  that  set  of  subsecs  of  B  which  contains  the  sub¬ 
set  M  of  B  If,  and  only  If,  there  exist  functions  F  and  G  In  D 

such  that  (F,*,G)  *  M  .  In  this  section,  we  want  to  do  two  things:  -  we 

would  like  to  show  how  the  theory  in  Chapter  II  of  ARL  65-75  suggests  a  way 

of  generating  special  subsets  of  the  non-empty  members  of  the  class  C£  by 
employing  the  group  of  all  one-one,  Lebesgue-meaoure-preserving  mappings  of 
the  closed  interval  I  ■  [0,1]  onto  itself;  -  and  we  would  like  to  pro¬ 

vide  illustrative  material  of  sufficient  complexity  to  justify  a  few  of  the 
elaborate  details  of  the  theory  which  was  constructed  around  class  in 

ARL  65-75.  In  order  to  realize  the  latter  intention,  we  shall  make  use  of 
some  pathological  functions  throughout  the  discussion;  and  we  shall  define 
these  functions  by  means  of  sprays.  The  subject  "Spray-forms  and  Sprays"  is 
considered  in  section  3  of  the  Appendix  in  ARL  65-75;  however,  the  brief 
discussion  of  this  subject  which  follows  should  be  adequate  for  our  present 
needs. 

A  sequence  (by  which  term  we  mean  an  infinite  sequence)  may  be  thought 
of  as  an  array  of  things,  called  its  terms,  which  have  been  entered  one 
after  another  into  the  places  of  a  figurative  structure  which  we  shall  call 
the  sequence-form.  If  the  entries  (i.e.,  the  terms)  of  a  particular  sequence 
are  real  numbers,  then  the  sequence  is  an  element  of  the  set  W  which  was 
defined  in  section  2.  However,  a  particular  sequence,  which  is  produced  by 
entering  into  each  place  of  the  sequence- form  a  sequence  of  real  numbers, 
one,  and  only  one,  of  which  is  non-zero  and  is  1  ,  is  a  matrix  in  the  set 
T  which  was  defined  in  section  2.  By  such  examples  as  the  latter,  one 
could  justify  the  idea  that  matrices  are  special  sequences.  For  this  reason, 
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we  do  not  hesitate  to  bypasa  matrices  in  our  search  for  a  generalization  of 
the  sequential  notion. 

For  each  natural  number  n  ,  a  term  in  the  fundamental  sequence  of 
natural  numbers,  let  ?(l,n)  denote  one,  and  only  one,  place  in  the 
sequence- form;  and,  for  each  place  in  the  sequence-form,  let  there  be  one, 
and  only  one,  natural  number  n  such  that  P(l,n)  denotes  that  place. 

Each  place  P(l,n)  in  the  sequence- form  has  one,  and  only  one,  successor : 
the  place  P(l,n+1)  ;  each  place  P(l,n)  which  is  different  from  P(l,l) 
has  one,  and  only  one,  predecessor :  the  place  P(l,n-1)  ;  and  the  place 
P(l,l)  has  no  predecessor. 

By  analogy,  we  produce  a  generalized  version  of  the  sequence-form 
which,  for  a  positive  integer  k  ,  is  called  the  k-fpoted  spray-form.  For 
each  ordered  pair  of  positive  integers  (m,n)  ,  in  which  mi  k"  *  ,  let 
P(m,n)  denote  one,  and  only  one,  place  in  the  k-footed  spray-form;  and, 
for  each  place  in  the  k-footed  spray-form,  let  there  be  one,  and  only  one, 
ordered  pair  of  positive  integers  (m,n)  ,  in  which  m  £  kn  *  ,  such  that 
P(m,n)  denotes  that  place.  Each  place  P(m,n)  in  the  k-footed  spray-form 
has  k  ,  and  only  k  ,  successors ,  the  t'th  of  which  is  the  place 
P(km-k+t,n+l)  ;  each  place  P(m,n)  which  is  different  from  P(l,l)  has 
one,  and  only  one,  predecessor ;  the  place  P(h,n-1)  ,  in  which  h  is  the 
greatest  integer  less  than  (m+k)/k  ;  and  the  place  P(l,l)  has  no 
predecessor. 

The  on»;-footed  spray-form  is  the  sequence  form.  The  two-footed  spray- 
form,  which  we  shall  call  simply  the  spray- form,  will  now  be  given  our 
attention.  Figure  1  is  an  illustrative  array  of  fifteen  of  the  symbols 
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which  denote  places  of  the  spray-form;  in  it,  the  fourteen  arrows  are 
directed  from  seven  of  the  place  symbols  toward  the  symbols  which  denote 
their  successors. 


Figure  1 
The  Spray-Form 

Pd.D 


P<3,3)  P(4,3) 

/\  A 

P(5,4)  P(6,4)  P(7,4)  P(8,4) 

We  shall  refer  to  the  first  and  second  successors  of  a  place  as  its  left 
and  right  feet,  respectively;  and  we  shall  refer  to  the  predecessor  of  a 
place  as  its  head.  The  place  P(m,n)  is  said  to  be  in  the  m' th  position 
on  the  n'th  level  of  the  spray-form.  The  place  P(l,l)  is  called  the 
topmost  place.  A  sequence  of  places  of  the  spray-foim,  in  which  the  place 
P(m,n)  is  the  first  term  and  each  term,  a  place,  is  followed  immediately 
by  either  its  left  foot  or  its  right  foot,  is  said  to  be  a  limb  of  P(m,n)  . 
A  limb,  all  of  whose  terms  after  the  first  term  are  left  (or  right)  feet, 

Is  called  a  left  limb  (or  a  right  limb) ;  and  a  limb  which  is  either  a  left 
limb  or  a  right  limb  is  called  a  straight  1 imb.  Thus,  eacl.  place  of  the 
spray-form  has  exactly  two  straight  limbs:  a  right  limb  and  a  left  limb. 

The  left  limb  of  a  place  P(m,n)  which  is  not  a  term  of  the  left  limb 
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of  P(l,l)  has  as  its  so-called  heart  that  place  whose  right  foot's 
left  limb  has  P(m,n)  as  one  of  its  terms;  the  right  limb  of  a  place 
P(m,n)  which  is  not  a  term  of  the  right  limb  of  P(l,l)  has  as  its  heart 
that  place  whose  left  foot's  right  limb  has  P(m,n)  as  one  of  its  terms; 
the  left  limb  of  a  place  which  is  a  term  of  the  left  limb  of  P(l,l)  has 
no  heart;  and  the  right  limb  of  a  place  which  is  a  term  of  the  right  limb 
of  P(l,l)  ha*  no  heart.  The  left  and  right  limbs  of  P(4,4)  have  the 
respective  hearts  P(2,3)  and  P(l,l)  ;  the  left  and  right  limbs  of  P(81,8) 
have  the  respective  hearts  P{3,3)  and  P(41,7)  ;  the  right  limb  of 
P(l,3)  has  the  heart  P(l,2)  ;  and  the  left  limb  of  P(l,3)  has  no  heart. 

By  the  branch  of  a  place  P(m,n)  is  meant  that  part  of  the  spray-form 
which  contains  a  place  P(a,b)  if,  and  only  if,  P(a,b)  is  a  term  of  a 
limb  of  P(m,n)  .  A  branch  is  structurally  the  same  as  the  spray-form. 

A  proper  limb  of  the  spray-form  is  any  limb  of  the  topmost  place  which 
is  not  the  left  limb  of  P(l,l)  ,  and  which  has  infinitely  many  terms  that 
are  left  feet. 

A  spray  (for  "spreading  array")  is  envisioned  as  being  the  array  which 
is  produced  by  entering  objects,  called  entries  „  into  the  places  of  the 
spray-form;  this  is  analogous  to  the  treatment  of  a  sequence  as  the  result 
of  having  entered  objects,  called  terms,  into  the  places  of  the  sequence- 
form. 

By  a  limb  L*  in  £  spray  S  (which  is  not  the  same  as  a  limb  of  a 
place  in  the  sprcy-form)  is  meant  the  sequence  of  entries  in  the  terms  of 
a  proper  limb  L  of  the  spray-form;  the  limb  L1  in  S  is  said  to  occupy 
the  proper  limb  L  of  the  spray-form. 

We  shall  have  need  to  consider  only  those  sprays  whose  entries  are 
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real  numbers,  and  whose  limbs  are  convergent  sequences.  Each  such  spray  S 
is  called  a  valuable  spray ;  and  the  limit  of  each  limb  in  S  is  said  to  be 
the  value  of  the  limb.  If  S  is  a  valuable  spray,  then  S  Is  said  to 
spread  over  that  subset  M  of  th.  real  line  R  which  contains  the  real 
number  x  if,  and  only  if,  there  is  a  limb  in  S  whose  value  is  x  •  If 

S  is  a  valuable  spray,  and  if  no  two  limbs  in  S  have  the  same  value, 

then  S  is  said  to  be  a  real  spray . 

J  If  S  is  a  real  spray  which  spreads  over  R  ,  and  if  T  is  a 

yaluable  spray,  then  the  symbol  S  — ►  T  denotes  that  function  which  maps 
R  into  R  in  such  a  way  that,  at  each  point  x  of  R  ,  its  value  y  in 
R  is  the  value  of  that  limb  in  T  which  occupies  the  same  proper  limb  of 
the  spray-form  that  is  occupied  by  the  limb  in  S  whose  value  is  x  . 

We  are  now  ready  to  define  five  special  valuable  sprays.  This  is  done 
for  the  purpose  of  facilitating  our  definition  of  some  pathological  Baire 
functions  in  B  . 

The  spray  X  is  constructed such  a  way  that  each  of  its  entries  is 
a  ratio  of  two  integers  whose  greatest  common  divisor  is  1  .  For  each 
| -iositive  integer  n  ,  its  entries  in  the  places  P(l,n)  and  P(2n  \n)  of 
the  spray-form  are  — —  and  ,  respectively.  The  sequence  of  numerators 

of  its  entries  in  the  places  of  a  straight  limb  of  the  spray-form  which  has  a 
heart  is  an  arithmetic  progression  whose  common  difference  is  the  numerator 
of  its  entry  in  that  heart;  and  the  sequence  of  denominators  of  its  entries 
in  the  places  of  a  straight  limb  which  has  a  heart  is  an  arithmetic 
progression  whose  common  difference  is  the  denominator  of  its  entry  in  that 
heart.  The  first  five  levels  of  spray  X  are  shown  in  Figure  2. 
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Figure  2 
Spray  X 
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The  spray  X  is  constructed  in  such  a  way  that,  for  any  positive  inte¬ 
gers  m  and  n  which  are  such  that  2n  2  <  a  i  2°  1  ,  its  entry  in  P(m,n) 
is  the  square  of  the  entry  of  spray  X  in  P(m,n)  ;  and,  for  any  integers 
m  and  n  which  are  such  that  1  <  m  <  2n  *  ,  its  entry  in  P(m,n)  is  the 
negative  of  its  entry  in  P(2n  ^+l—tn,n)  .  The  first  five  levels  of  spray  X 
are  shown  in  Figure  3. 


Figure  3 
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The  spray  Y  la  constructed  In  such  a  way  that  each  of  Its  entries 


Is  a  ratio  of  two  Integers  whose  greatest  common  divisor  Is  1.  Its 

entries  In  P(l,l)  ,  P(1.2)  ,  P(2,2)  and  P(2,3)  are  |  ^  , 

27  X  5 

22  ^  »  re*Pectlvely.  For  each  Integer  n  >  2  ,  Its  entry  In 

P(l,n)  Is  l/2n+^  ;  and  Its  entry  In  each  place  of  the  branch  of  P(2,n+1) 

la  l/2n+*  .  If  -  Is  Its  entry  In  any  place  of  the  branch  of  P(2,3)  , 


then  its  entries  In  the  left  and  right  feet  of  the  place  occupied  by 
m  4ia9  4*+3  g 

£  are  ■  and  »  respectively;  and,  if  —  is  Its  entry  In  any 

place  of  the  branch  of  P(2,2)  ,  then  its  entries  in  the  left  and  right 

,  _  t  a  4a-9-V2b  ,  4a+3+"\/2b 

feet  ot  the  place  occupied  by  —  are  - ^ -  and  - — -  , 

respectively.  The  first  five  levels  of  Y  are  shown  in  Figure  4. 


Pigure  4 
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The  spray  Z  is  constructed  in  such  a  way  that  its  entry  in  each 

place  of  the  branch  of  P<2,2)  is  a  ratio  of  two  integers  whose  greatest 

cocsaon  divisor  is  1;  and  its  entry  in  each  place  which  does  not  belong 
to  the  branch  of  P(2,2)  is  0.  Its  entry  In  P(2,2)  is  ^  ;  and,  if 

is  its  entry  in  any  place  of  the  branch  of  P(2,2)  ,  then  its  entries 

D 
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in  the  left  and  right  feet  of  this  place  are 
respectively.  The  first  five  levels  of  Z 


2a-l 


and 


2a+l 


2b  2b 

are  shown  in  Plgure  5. 


» 


Figure  5 
Spray  Z 
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The  spray  V  Is  constructed  In  such  a  way  that  Its  entry  In  each 
place  of  the  branch  of  p(l,2)  Is  0  ;  and  Its  entry  In  each  place  which 
does  not  belong  to  the  branch  of  P(l,2)  Is  a  ratio  of  two  Integers  whose 
greatest  common  divisor  is  1  .  Its  entries  In  P(l,l)  ,  P{2,2)  and 


P(3,3)  are  -  ,  —  and  *  respectively.  If  ^  is  its  entry  in  any 

place  of  the  branch  of  P(3,3)  ,  then  its  entries  in  the  left  and  right 

,  _  ,  ,  4a-27- V5b  ,  4a+9+V2b  .  , 

feet  of  this  place  are  - — -  and  - —  ,  respectively.  For  any 

integer  n  >  2  ,  the  entries  of  V  in  P(2nl,n)  and  P(2n-l,n+l)  are 

/02n+l  »n- 1  .  2n+l  ,  _2n+3  ,  ,n+l  ._2n+3  .  .  . 

(2  -3*2  - 1 ) / 2  and  (2  -5*2  -13)/2  ,  respectively;  and, 

i  n 

if  ^  is  its  entry  in  any  place  of  the  branch  of  P(2  -l,n+l)  ,  then  its 

4a-9-"V2b 

entries  in  the  left  and  right  feet  of  this  place  are  - — —  and 

__  4b 

4a+3+v2b 

- 4b - *  re8P*-t*-vely •  The  fir*t  five  levels  of  spray  V  are  shown  in 

Figure  6. 
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Figure  6 
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By  making  use  of  the  sprsys  X  ,  Y  ,  Z  ,  X  and  V  ,  the  Baire 
functions  £  ,  £  and  G  in  B  are  defined  as  follows:  £  -  X  — ►  Y  , 

£  •  X  — *>  Z  and  G  -  X  — >  V  .  Now  that  the  letters  X  ,  Y  ,  Z  and  V 
have  served  their  purpose  in  the  definition  of  the  special  functions  £  , 

£  and  G  ,  we  free  them  for  other  uses  in  this  discussion. 

Some  particular  values  of  these  special  functions  are:  F(-3/2)  ■  1/16  , 

* 

£(-101/100)  -  1/16  ,  £(-l )  -  3/16  ,  £(-2/3)  -  51/256  ,  £((l-V5)/2)  -  1/5  , 

1(0)  -  5/8  ,  £(V?)=  74/85  „  £(3/2)  -  451/512  ,  7(0)  -  0  ,  F(l)  -  1/8  , 

£(V2)  -  3/20  ,  £(3/2)  -  5/32  ,  G(-l/i00)  -  0  ,  G(0)  -  3/8  , 

G(4/25)  -  909/2043  ,  C(9/16)  -  1181/204 8  ,  G (1)  -  25/32  and 

G(3/2)  -  3439/4368  . 

Throughout  the  remainder  of  this  section,  we  shall  have  occasion  to 
repeat  some  of  the  definitions  given  previously.  This  will  be  done  in  the 
process  of  providing  the  kind  of  illustrative  material  which  we  hope  will 
support  the  procedure  that  was  followed  in  the  proof  of  the  theorem  in 
Chapter  II  of  ARL  65-75. 
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The  symbol  Ej  denotes  the  class  of  all  Borel  sets  In  R  .  The  symbol 
B  denotes  the  set  of  all  Balre  functions  which  map  R  into  R  .  The  symbol 
I  denotes  the  closed  Interval  tO,l]  in  R  .  The  symbol  D  denotes  the 
set  of  all  distribution  functions  which  map  R  into  I  .  For  example,  the 
special  functions  F  ,  F  and  G  are  in  B  ;  and  the  functions  F  and  G 
are  in  D  • 

For  any  h  €  B  ,  the  symbol  h  denotes  that  function  which  maps  R 
into  B^  in  such  a  way  that  its  value  at  each  y  €  R  is  the  Borel  set 
h(y)  ■  (x  :  h(x)  i  y)  .  The  special  symbol  W  denotes  that  function  in 
B  whose  value  at  each  x  (  R  is  M(x)  -  x  ;  therefore,  for  each  x  C  R  , 
to>(x)  denotes  the  half-line  (-**,xl  . 

For  any  F  €  D  ,  the  symbol  denotes  that  probability  measure 

which  maps  B1  into  I  in  such  a  way  that,  for  each  x  €  R  ,  the  P  - 

meaaure  of  the  Borel  set  5(x)  is  PjJLwCx)]  »  F(x)  .  For  any  F  £  D  and 

h  €  B  ,  there  exists  one,  and  only  one,  function  G  in  D  whose  value  at 

each  y  £  R  is  G(y)  *  PpLh(y)}  ;  and  the  fact  that  the  functions  F  ,  h 

and  G-  are  related  to  one  another  in  this  way  is  stated  symbolically  by 

"  (F,h,G)  For  any  F  €  D  and  G  €  D  ,  the  symbol  (F,*,G)  denotes  that 

subset  of  B  which  contains  h  €  B  if,  and  only  if,  (F,h,G)  is  a  true 

statement.  For  example,  the  special  set  (F,*,C)  contains  the  special 

2 

function  h  in  B  whose  value  at  each  x  £  R  is  h(x)  -  x  .We  intend 
to  describe  a  method  of  obtaining  various  functions  in  a  set  (F,*,  G) 
when  one  such  function  is  known. 

For  any  F  £  D  ,  the  symbol  R(l;«*)  denotes  the  set  of  all  points  of 
discontinuity  of  F  ;  and  the  symbol  R'(F,-«»)  denotes  the  difference  set 
R  -  R (F ;«*)  .  For  example,  R(F;«)  contains  x  €  R  if,  and  only  if,  x  is 
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either  e  negative  integer  or  a  rational  number  in  the  open  half-line  (-1  ,«•)  ; 
and  K(C;w)  contain*  y  (  R  if,  and  only  if,  y  la  the  aquare  of  a  non- 
negative  rational  number. 

For  any  F  €  D  ,  the  aymbol  (F;oo)  denotea  that  aubset  of  R'(F;s») 

which  contain*  x  6  R'(F;«a)  if,  and  only  if,  for  each  real  number  r  <  x  , 

F(r)  <  F(x)  ;  and  the  aymbol  R^(F;«»)  denote*  the  difference  aet 

l' (F;*»)  -  Rc(F;*»)  .  Furthermore,  for  each  x  €  R  ,  the  aymbol  Rc(F;x) 

denote*  the  interaection  aet  Rc(F;»)  A  e(x)  ;  and  the  aymbol  Rc(F;-e») 

denote*  the  empty  act  t  .  If  R^(F;«*)  4  <t  ,  then  it  la  either  a  aingle 

Interval  with  no  left  endpoint  or  the  union  of  pairvlae  disjoint  intervale 

with  no  left  endpoints;  and,  over  each  such  interval  in  the  composition  of 

R^(F;w)  ,  the  function  F  has  only  one  value  (l.e.,  it  has  a  so-called 

"plateau").  For  example,  R^(F;«)  ia  the  set  of  all  non-integers  of  the 

half-line  4>(-l)  .  Despite  the  fact  that  the  arbitrary  value  F(x)  of  a 

function  F  €  D  increases  as  x  increases  over  R  (V;°o)  ,  it  need  not  be 

c 

true  that,  if  X  is  a  Borel  subset  of  Rc(F;»)  whose  Lebesgue  measure  is 
positive,  then  Pp£X]  is  positive.  It  happens,  for  example,  that,  if  X  is 
that  subset  of  Rc(F;»)  which  consists  of  all  the  Irrational  numbers  in  the 
open  Interval  (-1,0)  ,  then  P^IXj  ■  0  . 

Let  R  denote  that  special  set  which  consists  entirely  of  all  the 

poluta  of  R  and  the  two  special  non-real  points  ao  and  -oo .  For  any 

F  €  D  ,  let  the  symbol  ?  denote  the  function  which  maps  R  into  I  in 

such  a  way  that  its  value  at  each  point  x  of  R  is  P?[Rc(F;x)]  .  For 

example,  if  F  ia  F  ,  then,  for  each  x  €  R  »  7(x)  ■  1F(x)  . 

For  any  F  €  D  ,  the  symbol  R^(F;«)  denotes  that  subset  of  Rc(F,vo) 
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which  contains  x  •  (F;»)  If,  and  only  If,  for  each  real  number  r  <  x  , 

'p(r) <  'f(x)  .  For  example,  R^(F;ao)  la  the  set  of  all  positive  lrrattooal 
numbers;  and  the  set  RW(G;«)  Is  tne  set  of  all  positive  real  numbers  which 
are  not  the  squares  of  rational  numbers. 

For  any  F  t  D  ,  let  the  symbol  denote  the  closed  Interval 

onto  which  the  set  R  Is  mapped  by  the  function  ?  ;  let  the 
symbol  F  denote  that  function  which  maps  RW(F;«)  Into  Ip  In  such  a 
way  thai”  Its  value  at  each  x  *  RW(F;«$  Is  F(x)  ■  F(x)  ;  let  the  symbol 

(  1^  If,  and  only  If, 

A  A 

of  F  Is  F(x)  •  u  ;  and 

let  the  symbol  J  denote  that  subset  of  T  which  Is  such  that 
F  F 

Jp  O  Hy  ■  Ij,  and  J^A  t  .  Since  F  Is  a  one-one  mapping  of  R^(F;«o) 

onto  Hp  ,  It  has  a  unique  Inverse  which  nay  be  denoted  by  the  symbol  F  *  . 
Since  the  Lebesgue  measure  of  Jp  Is  0  ,  the  Lebesgue  measure  of  Hp  la 
the  length  7(oo)  of  the  Interval  Ip  .  For  example,  since  Is  the 

closed  Interval  [0,1/4]  ,  and  since  is  that  subset  of  Ip  which 

contains  0  and  a  positive  real  number  u  if,  and  only  If,  there  exist 
Integers  m  and  n  which  are  such  that  0<m£2n^+^  and 
(2m-l)/2n+*  ■  u  ,  It  follows  that  Jp  contains  a  denumerably  Infinite 
number  of  points,  so  that  Its  Lebesgue  measure  Is  0  ;  consequently,  the 
Lebesgue  measure  of  both  the  set  Up  ■  Ip  —  Jy  and  the  Interval  Ip  la 
the  latter's  length  1/4  . 

If  a  is  a  one-one  mapping  of  a  set  S  onto  itself,  and  If  X  la  a 
subset  of  S  ,  then,  by  the  modified  restriction  of  <3  to  the  subset  X 
of  S  ,  we  mean  the  unique  one-one  mapping  g  of  X  onto  itself  which  Is 
defined  as  follows:  let  Y  -  S  -  X  ;  let  X^  be  that  subset  of  Y  which 
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Hp  denote  that  subset  of  Ip  which  contains  t 
there  exists  an  x  C  R  (F;oo)  at  which  the  valui 


contains  y  «  Y  If,  and  only  if,  a_i(y)  -  x  is  an  clement  of  X  ,  where 
O'1  denotes  the  unique  Inverse  of  a  ;  for  each  positive  Integer  n  ,  let 
X^  be  that  subset  of  X  which  contains  x  €  X  if,  and  only  if,  a(x) 
is  an  element  of  X^n  ;  let  X1  be  chat  subset  of  X  which  contains 

x  «  X  if,  and  only  if,  there  exists  a  positive  integer  n  such  that 

x  €  X^  ;  let  X"  ■  X  —  x’  ;  and  let  &  be  that  function  over  X  whose 

value  at  each  x  €  X  la  either  P(x)  -  x  or  g(x)  -  Q(x)  according  as 
x  f  X1  or  x  €  X1'  .  Since,  for  each  non-negative  integer  n  ,  and  for 

each  x  €  X^  ,  Of  *(x)  is  either  an  element  of  Y  or  an  element  of  the 

subset  X(n+1)  of  X  ,  there  can  be  no  more  elements  in  than  there 

are  in  X^  ;  consequently,  if  S  is  a  Borel  subset  of  R  ,  if  Ot  is  a 
one-one,  Lebesgue-measure-preservlng  mapping  of  S  onto  S  ,  if  X  is  a 
Borel  subset  of  S  ,  and  if  Y  -  S  -  X  is  a  set  with  Lebesgue  measure  0  , 
then  every  Borel  subset  of  X*  has  Lebesgue  measure  0  . 

Let  M  denote  the  set  of  all  one-one,  Lebesgue-measure-preserving 
mappings  of  the  closed  Interval  I  onto  itself.  For  any  a  €  M  and  (3  6  M  , 
let  the  symbol  0a  denote  that  function  in  M  whose  value  at  each  u  €  I 
is  Pa(u)  »  P(Cr(u))  ;  and  let  the  symbol  a  *  denote  the  inverse  of  a  . 

The  set  M  constitutes  a  group  with  respect  to  the  binary  operation  that  is 
indicated  when  two  symbols  which  denote  functions  in  M  are  placed  side  by 
side  in  order  to  form  a  composite  symbol  which  denotes  a  function  in  M  . 

For  any  a  (  M  and  F  €  D  ,  let  the  symbol  Q denote  that  one-one, 

Lebesgue-measure-preserving  mapping  of  tue  closed  interval  1^  onto  i^, 

whose  value  at  each  u  t  I  is  either  Q  (u)  -  u  or  a  (u)  ■  ?(oo)a(u /?$»)) 

r  r  F 

according  as  ?(oo)  is  or  is  not  0  ;  and  let  the  symbol 
modified  restriction  of  to  the  subset  of  I  . 


denote  the 


'} 


f 


Lee  T  denote  that  function  which  maps  M*D*B  into  B  in  such  a 

way  that,  for  each  ordered  triple  (a,F,h)  in  MxD*B  ,  the  image  in  B 

of  (Ct,F,h)  is  the  function  TgFrh  whose  value  at  each  x  €  R  is  either 
TaF:h(x)  "  ^(F*1  (Op(F(x))))  or  TCfF:h(x)  »  h(x)  according  as  x  is  or  is 
not  a  point  of  Rw(F;oo)  .  The  image  in  B  of  (a,F,h)  C  M*D*B  under  the 

mapping  T  is  denoted  by  the  symbol  T^F : h  in  order  that  one  may  be  able 

to  denote  some  other  functions  by  convenient  symbols.  Thus,  for  any  Or  «  M  , 

the  symbol  TQ  denotes  that  function  which  maps  Ox  B  into  B  in  such  a 

way  that,  under  TQ  ,  the  image  in  B  of  the  ordered  pair  (F,h)  6  D«B 

is  T^Frh  ;  and,  for  any  ordered  pair  (a,F)  in  Mx  D  ,  the  symbol  (or 

operator)  TQF:  denotes  that  function  which  maps  B  into  B  in  such  a 
way  that,  under  T^F:  ,  the  image  in  B  of  the  function  h  6  B  is  TQF:h  . 

If  (Or,F,h)  is  any  ordered  triple  in  MxDxB  ,  and  if  G  is  that 

function  in  D  which  makes  (F,h,G)  a  true  statement,  then  the  subset 
(F,*,G)  of  B  contains  both  h  and  T^Frh  . 

For  each  ordered  triple  (a,b,c)  in  1*1x1  which  is  such  that 
a  £  b  and  b-a  ifi  c  ,  let  the  symbol  (a:b:c)  denote  that  function  in  M 

whose  value  at  each  u  €  i  is  either  (a:b:c)(u)  »  a+c-u  or  (a:b:c)(u)  »  u 

according  as  u  i»  or  is  not  a  point  of  the  union  of  closed  Intervals 
[a,b]  U  [ a-b+c,c3  . 

For  our  illustration  of  this  method  of  using  an  operator  T^:  in 
order  to  obtain  a  function  in  (F,*,G)  which  differs  from  h  €  (£,*,£) 


over  a  subset  of 

R  with  positive 

P^— measure,  we  choose  a 

to 

be  the 

convenient  mapping 

(1/4 : 1 / 3 :3/4) 

in  M  .  Then,  since  h  is 

that 

function 

in  B  whose  value 

at  each  x  f  R 

2 

is  h(x)  -  x  ,  the  function 

k 

-  TgFrh 
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1 


2 

la  B  it  such  that  ite  value  at  each  x  6  R  la  either  k(x)  -  1/x 
2 

k(x)  "  x  according  aa  x  la  or  la  not  an  irrational  number  in  the 
of  cloaed  Intervale  £l/2,  (V5-l)/2]  U  [(V5+l)/2,2]  . 


or 

union 
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III.  ROBUSTNESS  OP  TESTS 


1.  An  Illustration  of  the  Notion  of  Robustness  of  a  Test 

In  order  to  put  a  particular  notion  of  robustness  of  tests  into  a 
somewhat  realistic  and  tractable,  hypothetical,  experimental  situation, 
we  present  the  following  illustration: 

The  officer  in  charge  of  a  military  radar  operators'  school  wishes 
to  consider  a  change  in  the  training  regimen  that  has  been  followed  for 
several  years  to  one  which  is  much  less  costly.  A  group  of  n  entrants 
is  put  into  an  experimental  program  for  the  full  duration  of  the  course 
of  training.  It  has  been  the  practice  at  this  school  to  assign  to  each 
trainee  a  proficiency  rating  at  the  completion  of  his  training.  A 
recent  study  of  the  school's  records  seems  to  justify  the  assumption 
that  the  proficiency  ratings  are  normally  distributed.  Therefore,  next 
to  each  such  rating  in  the  school's  records,  but  not  in  the  service 
records  of  the  graduates,  is  placed  a  standardized  rating,  which  is  that 
transformation  of  the  given  proficiency  rating  that  adjusts  the  aggregate 
of  all  ratings  to  conform  with  the  assumption  that  they  are  normally  dis¬ 
tributed  with  mean  0  and  variance  1  .  It  is  Intended  that  the 
proficiency  ratings  of  the  trainees  in  the  experimental  group  shall  be 
transformed  (l.e.,  shall  be  standardized)  by  the  same  function  that  is 
used  on  all  the  proficiency  ratings.  And,  in  order  to  meet  the  require¬ 
ments  of  the  Analysis  of  Variance  methodology  in  testing  the  worth  of 
the  proposed  change  in  the  course,  it  is  assumed  that,  if  for  each  posi¬ 
tive  Integer  1  <  n  ,  is  the  random  variable  whose  values  are  the 
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possible  standardized  ratings  that  the  l'th  trainee  In  the  experimental 
course  might  achieve  upon  the  completion  of  his  training,  then  the  vari¬ 


ance  of  X^  is  1  .  Two  further  reasonable  assumptions  are  made:  if 
i  jl  J  ,  then  Xt  and  Xj  are  independent  random  variables;  and,  for 
each  positive  integer  i  <  n  ,  X^  is  normally  distributed  with  mean  p 
The  distribution  function  and  the  probability  density  function  of  X^  , 
which  la  N(p,l)  ,  are  denoted  simply  by  and  ^  ,  respectively; 

and  and  ^  may  be  denoted  simply  by  J  and  ^  ,  respectively. 

Furthermore,  the  probability  measure  which  induces  i  i«  denoted  by  P  • 
Thus,  the  sample  mean  of  the  standardized  ratings  of  the  experimental 


Xj  +  ...  +  X 

units  is  the  normal  random  variable  X  -  - - - -  with  mean  p  and 

variance  -  . 

n 

Now,  a  real  number  <2  in  the  Interval  (0,1)  is  chosen;  and,  for 
each  real  number  v  ,  the  symbol  c(v,0()  denotes  that  real  number  which 
satisfies  the  following  equation:  p{X  <  c(v,a)  ;  p  -  v)  ■  a  .  The  num¬ 
ber  c(v,a)  is  the  right  endpoint  of  an  open  half-line  critical  region 

r  . 


With  all  these  things  under  consideration,  the  question  of  whether 
or  not  the  regular  training  program  should  be  replaced  by  the  experimental 
training  program  is  taken  up  in  the  folloving  test  which  te3ts  :  p  ■  0 

against  the  :  p  <  0  . 

If  X  <  c(0,a)  ,  then  reject  the  null  hypothesis  Hq  ; 

and  if  X  >  c(0,a)  ,  then  do  not  reject  Hq  . 

The  power  function  of  this  test  is  defined  at  each  real  point  v  by 

P(v)  -  P{X  <  c(0,a)  :  p  -  v}  -  P  <  -c(2>a)~v  ]  -  j  (s.(9fa)-v  \  . 

U/N/n  1/Vn  J  '  1/Vn  ■ 
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Since  a  -  p(x  <  C(0,a)  :  p  -  0)  -  p(x  Vn  <  c(0,a)  Vn}  -  f  (c(0,O)  \/n)  , 

the  reel  nunsber  T  -  c(0,a)  is  f  ~1<p)/\Jn  ,  so  that,  at  each  real 
point  v  ,  the  value  of  (3  is  0(v)  ■  Jc  J  *  (a)  *  v\/n) 

Now,  euppose  that  the  underlying  assumption  that  the  variance  of 
X^  is  1  is  changed;  that  is,  suppose  that,  for  each  positive  integer 
i  <  n  ,  X^  is  N(n,cr  )  •  Under  this  assumption,  X  is  N(^,cr  /n)  . 

Then  too,  the  power  function  of  the  above  test,  in  which  the  critical 

region  P  is  the  open  half-line  with  right  endpoint  T  ,  has  as  its 
value  at  each  real  point  v  the  number 


Vv> 


p{x  <  r  :  M  -  v)  -  p-^ 
V  a/  Vn  '  V 


X-v  <  Jf 


<r/ Vn 
-1 


cr/  Vn 
(Ot)-vVn  ^ 

a  )  ' 


For  convenience,  the  functions  5  and  A  are  defined  at  each  real 
point  v  at  follows:  6(v)  ■  ^(v)  -  P(v)  and  A(v)  -  |k(v)|  •  The 
first  derivative  of  6  la  6'  whose  value  at  each  real  point  v  is 

8'(v)  -  -  -^2  ■  )  +  Vn  (^(  |  _1(Q)  -  v  Vn  )  . 


For  cr  /  1  ,  6*  (v)  -  0  if,  and  only  if,  v  is  either 


— 

1  1  " 

1 

t  -1 _ „  1 2  In  cr 

1  (“>  *  CT,  2 

1 

or  — — 

1  -1 .  1  2  In  o 

3  (“)  +  2 

Vn 

V  *  -i 

Vn 

V  <TZ-1 

Therefore,  for  each  real  number  v 

Tin 


0  <  A(v)  < 


I 
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where  the  upper  bound  on  the  right  ia  useful  only  if  cr  la  close  to  1  . 
In  a  rough  way.  It  can  be  said  that,  with  respect  to  the  above  mentioned 
change  in  the  underlying  assumptions,  the  measure  of  robustness  of  the 
test  under  discussion  is 


i 


This  measure  is  greater  than  1/2  (i.e.,  the  test  is  more  than  50Z 

robust)  if  c r  is  in  the  interval  (e  \e)  . 

i  i 

I 

2.  A  Generalization  of  the  Notion  of  Robustness  of  a  Teat 

In  the  last  section,  we  tried  to  illustrate  a  way  of  measuring  the 
robustness  of  a  test.  By  considering  a  contrived  experimental  situation 
and  a  convenient  test,  we  produced  a  number  between  0  and  1  which  was 
to  serve  as  a  rough  indicator  of  the  ineffectiveness  of  a  particular  change 
in  the  assumptions  underlying  the  test  as  a  disturbing  influence  on  its 
power  function.  This  number,  which  could  be  expressed  as  a  percentage, 
was  called  the  measure  of  robustness  of  the  test  with  respect  to  the 
particular  change  in  its  underlying  assumptions. 

I 

In  this  section,  we  shall  give  a  skeletal  version  of  a  generaliza¬ 
tion  of  the  illustration  in  the  last  section.  In  doing  this,  we  begin 
with  a  review  of  some  notational  conventions  which,  in  somewhat  more 
detail.,  are  treated  in  ARL  65-75. 

Por  each  positive  integer  n  ,  let  R^  denote  the  Euclidean  space 
of  n  dimensions;  for  each  point  x  of  Rr  and  each  positive  integer 
t  <  n  ,  let  x£  denote  the  t’th  coordinate  of  x  ;  let  R  -  R^  ; 
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let  B  denote  the  class  of  all  Borel  subsets  of  R  ;  for  each  point 
n  n 

x  of  R  ,  let  ( denote  that  member  of  B  which  contains  the 
n  n 

point  u  of  R^  if,  and  only  If,  for  each  positive  Integer  C  5  n  » 

u  <  x  ;  let  D  denote  the  set  of  all  distribution  functions  which 
t  -  t  n 

map  R  into  the  closed  interval  I  -  fO,l]  in  R  :  let  M  denote 
n  n 

the  set  of  all  probability  measures  which  map  Br  into  I  ;  for  each 

F  £  Dn  ,  let  Pr(F)  denote  that  probability  measure  P  in  Mn  which 

induces  F  ,  so  that,  for  each  point  x  of  Rr  ,  p{u  :  u  €  [-«,x])  -  F(x) 

and,  for  each  P  €  ,  let  Df(P)  denote  that  distribution  function 

F  in  Dn  which  is  induced  by  P  ,  so  that  Pr(F)  ■  P  .  Furthermore, 

for  each  ordered  pair  of  positive  integers  (m,n)  ,  let  B  denote  the 

*is  n 

set  of  all  Baire  functions  which  map  R  into  R  ;  and,  for  each  ordered 

n  oa 

pair  (F,h)  in  the  Cartesian  product  sat  Dn  X  mBn  ,  let  (F,h,*)  denote 

that  subset  of  D  which  contains  the  function  G  in  D  if,  and  only 

m  m 

if,  for  each  point  y  of  R  ,  G(y)  *  P(x  :  h(x)  €  f-~,y]}  ,  where 

CD 

P  -  Pr (F)  . 


If,  for  positive  integers  1  and  j  ,  V  denotes  a  non-empty 
subset  of  R^  and  W  denotes  a  non-empty  subset  of  R^  ,  then  the 
symbol  V:W  denotes  that  subset  of  R^+^  which  contains  the  point  u 
of  *nd  °nly  if»  there  exist  points  v  of  V  and  w  of  W 

such  that,  for  each  positive  integer  t  <  i  +  J  ,  ufc  is  either  wt_£ 
or  v£  according  as  t  does  or  does  not  exceed  i  ;  and,  if  this  point 
u  is  in  V:W  because  such  points  v  of  V  and  w  of  W  do  exist, 
then  u  may  be  denoted  alternatively  by  the  symbol  v:w  .  If  VCR 


and  DcR  ,  then  V:W  is  the  Cartesian  product  set  V  x  W  ;  however. 


if  V  c  R  and  W  c  R^  ,  then  V;W  need  not  be  a  Cartesian  product 
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subset  of  R^  . 

If  S  is  soy  non-empty  set,  then  s  function  P  which  maps  the 
Cartesian  product  set  S  X  S  into  R  is  said  to  be  a  metric  for  S  if, 
and  only  if,  the  following  three  conditions  are  satisfied: 

(1)  if,  for  any  ordered  pair  (a,b)  of  S  X  S  ,  P(a,b)  -  0  , 

then  a  -  b  ; 

(2)  for  each  a  e  S  ,  P(a,a)  -  0  ;  and 

(3)  for  each  ordered  triple  (a,b,c)  of  S  *  S  X  S  , 

P (a, c)  +  P (b , c )  >  P(a,b)  . 

Consequently,  for  each  (a,b)  €  S  x  S  ,  P(a,b)  -  [p(a,b)  +  P(a,b)]/2  > 

P(a,a)/2  ,  which,  by  (2),  is  zero,  so  that  P(a,b)  is  non-negative; 

and,  furthermore,  P(a,b)  »  P(o,b)  +  P(a,b)  by  (2),  P(b,b)  +  P(a,b)  > 

P(b,a)  by  (3),  P(b,a)  »  p(a,a)  +  P(b,a)  by  (2)  and  P(a,a)  +  P(b,a)  > 

P(a,b)  by  (3),  so  that  P(a,b)  >  P(b,a)  >  P(a,b)  and  P(a,b)  »  P(b,a) 

If  the  range  of  values  of  P  is  a  subset  of  I  ■  [0,l]  in  R  ,  then 

P  is  said  to  be  a  metric  for  S  limited  to  I 

Now,  let  i  ,  j  ,  m  ,  and  u  be  positive  integers;  let  the 

non-empty  subsets  V  and  W  of  and  R ^  be  called  parameter 

spaces ;  let  the  member  T  of  be  called  a  critical  region;  and  let 

the  function  h  of  B  be  called  a  transformation.  For  each  point 

v:w  of  V:W  ,  let  F  be  a  unique  distribution  function  in  D  ; 

v:w  n  * 

let  P  -  Pr(F  )  ;  and  let  G  be  that  distribution  function  in 

v : w  v  iw  v  jv 

the  subset  (F  ,h,*)  of  D  .  For  each  point  w  of  W  ,  let  3  , 

'  v:w  m  r  w 

called  a  power  function,  be  that  function  which  maps  V  into  the  interval 

I  whose  value  at  each  point  v  of  V  is  Pw(v)  ”  '•  *UX)  e 

Let  S  be  that  set  of  functions  which  map  V  into  I  ,  which  contains 
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Che  function  f  If,  and  only  if,  there  exists  a  w  €  W  such  that 
P  »  f  .  Finally,  let  P  be  a  Bietric  for  S  United  to  I  .  For 

V 

example,  P(f»g)  could  be  the  least  upper  bound  of  |f(v)  -  g(v)|  , 

where  f  e  S  ,  g  £  S  and  v  €  V  . 

Suppose  that,  in  an  experimental  situation,  a  decision  is  to  be 
based  on  the  way  a  random  variable  X  whose  range  of  values  is  R 

n 

is  distributed.  It  is  thought  that  the  distribution  function  of  X  is 
Fr-a  >  where  r:s  e  V:W  .  A  test  of  the  hypothesis  that  r  Is  the 
proper  parameter  in  V  is  designed;  end  it  makes  use  of  the  transfor¬ 
mation  h  as  well  as  the  critical  region  T  which  is  such  that 
Pg(r)  is  some  small  number  C*  in  I  .  j,_.is  test  may  be  stated  as 
follows:  if  h(X)  e  T  ,  reject  ;  otherwise,  do  not  reject  Hq  . 

The  power  function  of  this  test  is  0^  ,  where  the  parameter  s  is  a 

fixed  point  of  W  and,  hence,  is  an  underlying  assumption  of  the  test. 

In  ascribing  to  this  test  some  measure  of  robustness  with  respect  to  a 
change  in  the  underlying  assumption  Concerning  the  parameter  s  , 
another  test,  with  h  and  T  unchanged,  but  with  s  changed  to  w  e  W 
yields  the  power  function  P  .  And  the  number  1  -  P(p  ,0  )  serves 

as  a  measure  of  the  robustness  of  the  test  with  respect  to  the  particular 
change  in  its  underlying  assumptions. 

Of  course,  the  utility  of  such  a  measure 'of  robustness  of  a  test 
will  depend  on  the  choice  and  general  acceptance  of  the  metric  P  . 

Even  though  these  measures  of  robustness  are  in  I  and  are  intuitively 
appealing,  it  is  not  easy  to  defend  the  choice  of  any  one  of  them  with 
general  considerations. 
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3.  The  Kolmogorov  Metric  or  the  Space  of  Location  Families  of 
Distribution  Function* 

In  order  to  define  the  notions  of  robustness  of  a  statistical 
test  and  robustness  of  a  transformation,  some  quantitative  measure  of 
the  distance  between  power  fui  ctlons  and  the  distance  between  distribu¬ 
tion  functions  appears  to  be  highly  desirable.  It  would  appear  that 
several  authors  (see,  e.g.,  [l]-[iO])  have  held  this  view  to  a  greater 
or  lesser  extent.  In  section  1,  an  example  was  given  which  exploits 
the  Kolmogorov  metric  on  the  space  of  normal  distribution  functions  as 
a  measure  of  the  robustness  of  a  test.  The  power  functions  being  dealt 
with  In  this  example  are  very  simply  related  to  distribution  functions, 
a  situation  which  will  not  occur  In  general.  However,  for  those  situa¬ 
tions  for  which  It  does  occur.  It  appears  that  the  assumption  of  normal¬ 
ity  does  not  play  a  large  part  In  determining  the  outcome.  The  present 
section  Is  an  attempt  to  Investigate  the  Kolmogorov  metric  on  the  space 
of  distribution  functions  of  location-parameter  families. 

We  begin  with  definitions  and  notation.  Let  I  be  a  subset 

of  the  real  line,  and  let  (F(x;n)  :  u  €  1}  be  a  family  of  distribution 
functions  on  the  real  line.  We  may  suppose  that  the  distribution  func¬ 
tions  are  right-continuous,  and  that  the  ordering  on  I  Is  that  Induced 
by  the  natural  order  on  the  real  line.  We  will  call  the  family 
(F(x;w)  :  W  e  I)  a  location-parameter  family  with  respect  to  I  if  and 
only  If  there  Is  a  distribution  function  G  such  that  for  each  real  x  , 
F(x;p)  -  G(x-p)  ,  for  each  u  €  I  •  It  is  obvious  that  euch  a  family 
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is  stochastically  increasing;  i.e.,  if  X  €  I  and  p  «  I  with  X  >  ^  , 

then  F(x;X)  <  p(x;p)  ,  tor  all  real  x 

Since  we  will  require  a  measure  of  distance  between  distribution 

functions,  we  Introduce  one  such  measure,  often  railed  the  Kolmogorov 

metric.  For  each  pair  p  ,  X  with  ^  e  I  and  X  e  I  ,  and  for  each 

pair  F  ,  F\  of  distribution  functions,  wt  set  d(x;p,X)  -  |f(x;u)  -  F(x;X)| 
M  a 

and  p (F  ,F^)  -  sup  J(x;p,X)  It  is  readily  verified  ;^at  p  is 

^  x 

ll  I 

Indeed  a  metric  on  the  space  of  distribution  functions.  As  noted  above, 
a  location-parameter  family  (with  respect  to  I  )  is  stochastically 
ordered.  In  view  of  this,  we  may  dispense  with  the  absolute  value  signs 
in  the  definition  of  d(x;p,X)  if  we  assume  that  X  >  ^ 

Because  of  the  exploratory  nature  of  our  study,  we  will  make  a  series 
of  Increasingly  stringent  assumptions  oj  the  location-parameter  family 
{F(x;p)  ;  p  €  1}  ,  and  at  each  stage,  study  the  effect  of  each  additional 

assumption  on  the  Kolmogorov  metric. 

Assumption  l.  G  is  a  continuous  function. 

^  Cono-  "lence.  d(x;u,)j  is  a  continuous  function  f  x  . 

Proof.  d(x;p,X)  »  F(x;p)  -  F(x;X)  -  G(x-p)  -  G(x-X)  . 

Assumption  2.  G  is  a  strictly  increasing  function. 

Consequence.  To  eac,h  positive  real  number  q  ,  there  correspond 
real  numbers  a  (’I)  and  b(q)  ,  with  a(q)  <  b  (H)  ,  such  that  x  <  a(*l) 
implies  d(x;p,X)  <  q  and  b{Tl)  <  x  implies  d(x;p,X)  <  q  . 

Proof.  From  X  >  jj  and  G  strictly  Increasing  it  follows  thst 
F(x;p)  ■  G(x-p)  >G(x-X)  »  F(x;X)  ,  for  all  real  x  .  The  remainder 
of  the  proof  follows  from  the  fact  that  G  is  a  distribution  function. 
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We  may  lunxcrlze  Che  consequences  of  assumptions  1  and  2  aa  follows: 
Remark.  If  G  is  a  continuous  strictly  increasing  function, 
then  d(*;n,X)  la  a  continuous  function  which  attains  its 
(absolute)  maximum  on  a  finite  Interval. 

Proof.  d(*;u,X)  is,  when  restricted  to  ta^-,),b(T)]  ,  a  continuous 

function  on  a  closed,  bounded  interval,  and  hence  it  attains  its  (absolute) 
maximum  on  . 

Assumption  3.  G(x)  -  1  -  G(-x)  ,  for  all  real  x  ;  l.e.,  G 

is  a  symmetric  distribution  function. 

Consequence.  d(*;u*M  is  symmetric  about  the  point  ^(y+X)  . 

Proof.  d(~(^+X)  -  x;^,X)  -  G(|(ji+X)-x-u)  -  G(|(b+X)-x-X) 

-  G(i(X-u)-x)  -  G(|(n-X)-x)  -  1  -  G(x  -  |(X-b))  +  C-(x  -  |(n-X))  -  1 

-  G(x  +  |(X-p))  -  G(x  +  i(n-X))  -  G(|(u+X)  +  X-m)  -  G(|(n+X)+  x-X) 

■  d(^(p+X)  +  x;n,X)  ,  for  every  real  x  . 

Assumption  4.  G  is  absolutely  continuous,  with  probability  density 
function  g 

Consequence.  d'(»;p,X)  satisfies  d'  +x;^,X)  - 

-U'(^(u+a)  -  x;p,A)  ,  for  all  real  x  . 

Proof.  For  each  real  number  x  ,  define  k(x;u,X)  by  k(x;p,X)  ■ 
d(^(u+X)  +  x;^,X)  .  From  the  consequence  of  assumption  3,  it  follows 

that  k(«;u,X)  is  an  even  function.  Since  the  derivative  (when  it 
exists)  of  an  ever  function  is  an  odd  function,  the  proof  is  complete. 

We  note  that  assumptions  3  and  4  together  yield  the  following  result. 
Remark,  g  is  an  even  function. 

Proof.  The  existence  of  g  follows  from  assumption  4  and  its 
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*  \ -'ikovs  from  assumption  3. 

Assumption  5.  g  Is  increasing  for  x  <  0  . 

Consequence,  g  is  decreasing  for  x  >  0  ,  and  hence  g  has  a 

single  u^ximum  at  x  -  0  . 

Proof.  This  is  a  simple  consequence  of  the  assumption  and  the  fact 
that  g  is  an  even  function. 

Assumption  6.  g  has  a  derivative  g'  ;  thus,  g'  =  G1'  . 

Consequence.  If  x  >  0  ,  then  g'  (x)  <  0  . 

Proof.  Since  g  is  decreasing  for  x  >  0  ,  it  follows  that  its 

derivative,  g'  ,  must  be  negative  for  x  >  0  . 

We  c«n  aow  summarize  the  consequences  of  these  assumptions  in  the  following 

way. 

Theorem.  If  G  is  an  absolutely  continuous  function  which  is 
symmetric  about  zero,  and  if  g"  exists  and  is  an  increasing 
function  for  x  <  0  ,  then  P(F^,F^)  "  2G(^(X-y))  -  1  ,  for 

X  >  n  . 

Proof.  Because  X  ,  P (T  ,F,)  ■  sup  [F(x;u)  -  F(x;X)]  .  It 

'  W  K  x 

is  clear  that  d c|(u+X)  ;u,X)  -  2G(j(X-vi))  -  1  ,  and  that  sup  Tf(x;u)  -  F(x;X)] 

■  d(^(p+X);u,X)  follows  from  these  facts:  (1)  d '  (-|  (p+X)  ;u,X)  - 

g(|(X-M))  -  g(|(p-X))  -0  .and  (2)  d”  (|(u+X)  ;u,X)  -  g'(|(X-n))  -  g'(|(u-X)) 

-  2g'(|(X-u))  <  0  . 

As  we  will  show  presently,  this  theorem  furnishes  sufficient,  but 
not  necessary,  conditions  under  which  the  distance,  as  messured  by  the 
Kolmogorov  metric,  between  two  members  of  s  particular  family  of  distri¬ 
bution  functions  can  be  computed  as  a  function  of  the  distance  between 
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their  Indices.  We  list  here  some  cases  of  Interest  that  are  Included 


In  the  above  analysis. 


1  x 


1.  normal;  I  ■  8^  ;  F(x;p)  -  (2*)  2  J  exp  t-  |(t-p)2]dt  . 

x 

2.  Double  exponential;  I  -  8^  ;  F(x;y)  -  |  J  exp  f-|t-p|]dt  . 


3.  Cauchy;  I  ■  8^  ;  F(x;p)  - 


4.  Student's  t  ;  1  ■  8^  ; 


1  2 
*  arctan  [l+(x-p)  ]  ,  x  <  p 

^  ^  arctan  [l+(x-p)2]  ,  p  <  x 


/  \  r  /  \  -l-i  x  n+1 

-  r( ^  )  [r(  f  )V5*J  f  [i4o“l(t-w)2]'  2  dt  , 


for  each  natural  number  n  . 


5.  Gamma;  I  ■  Bj  ;  Fa>p(x;p)  -  P*°[ T(a) ] "l J  (t-p)0*1  exp  [-P'1  (t-p))dt 

M 

for  fixed  real  positive  Ot  and  fixed  real  positive  P  ,  x  >  u 

We  turn  now  to  a  location-parameter  family  of  distribution  functions 
with  nonsymmetrlc  members.  Take  I  -  ,  and  for  each  pci  ,  let 


F(x;p)  - 


0  .  x  <  p 

1  -  e:tp(p-x)  ,  p  <  x 


It  is  clear  that  the  family  (F(x;p)  :  p  €  i)  is  a  location-parameter 
family  with  respect  to  I  .  However,  the  distribution  function  G 
for  which  F(x;p)  «  G(x-p)  for  each  real  x  is 


G(y)  - 


o  ,  y  <o 

1  -  exp(-y)  ,  0  <  y 
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which  la  not  aymnecric  with  reapect  to  any  real  number.  Nevertheless, 
we  my  still  prove  the  following  result. 

Theorem.  If  F(x;n)  and  G(y)  are  given  by  (1)  and  (2) 
respectively,  then  for  any  4  €  I  and  X  e  I  ,  with 
X  >  |j  ,  we  have  P(F  -  G(X-(i)  . 

Proof.  Ve  distinguish  three  distinct  and  exhaustive  cases: 

(1)  -  <  x  <  n  ,  (11)  v*  <  x  <  X  ,  a  .id  (ill)  X  <  x  <  -  . 

(1)  For  x  <  u  ,  ?(x;vO  -  F(x;X)  -  0  ;  thus  d(x;u,X)  -  0  , 
for  x  <  vi  . 

(11)  For  vi  5  x  5  *  «  d(x;4,X)  can  he  written  as  d(x;4,X)  - 
F(x;p)  -  0  -  1  -  exp(vi-x)  ,  which  exhibits  the  fact  that  for  4  <  x  <  X  , 
d(x;|i,X)  Is  a  strictly  monotone  Increasing,  continuous  and  bounded 
function  of  x  .  Hence  d(x;vi,X)  attains  Its  maximum  value  (which  Is 
clearly  G(X-p))  at  x  «  X  . 

(Ill)  We  complete  the  proof  by  considering  this  case,  since  for 
X  <  x  ,  d(x;v*,X)  m  e  X(eX-eti)  <  e  X(eX-e^)  »  1  -  exp(X-4)  -  G(X~vi) 

Let  us  atcempt  to  general ire  this  thf-eo  to  Include  location- 
parameter  families  whose  members  are  not  necessarily  symmetric.  Thus, 
we  suppose  that  for  some  subset  I  of  ,  we  have  a  family 

(F(x;vO  :  m  €  l)  of  distribution  functions  such  that  for  each  4  €  I  » 
F(x;u)  •  0  ,  for  x  <  4  .  If  we  further  suppose  that  the  family 

(F(x;4)  :  4  €  i)  Is  a  location-parameter  family  with  respect  to  I  , 
with  F(x;4)  -  G(x~4)  for  all  real  x  ,  then, as  before,  this  family  Is 
stochastically  Increasing.  But  we  do  not  assume  now  that  G  is  synmietrlc 
with  respect  to  any  real  number.  Without  adding  any  restriction,  we 
suppose  that  X  >  4  ,  and  we  observe  that  for  -•  <  x  <  4  ,  we  have 
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F (x;X)  -  F(x;n)  "0  ,  so  that  for  x  <  n  ,  d(x;n,X)  ■  0  . 

If  we  now  suppose  that  G  Is  a  continuous  function  which  Is 

strictly  Increasing  whenever  It  is  positive,  then  as  above,  corresponding 

to  each  positive  real  number  q  ,  there  Is  a  real  number  a(q)  such 

that  a(n)  <  x  Implies  d(x*,n,X)  <  q  .  When  d(x;p,X)  Is  restricted 

to  the  closed  interval  t^,a(q)l  ,  It  is  a  continuous  function  on  a 

closed  interval  (which  Is  also  bounded),  and  hence  it  attains  its 

(absolute)  maximum  on  [p,a(q)]  .  From  the  nature  of  d(x;iu,X)  for 

x  e  [p,X]  ,  we  may  conclude  that  sup  d(x;p,X)  -  sup  G(x-p)  ■  G(X-p) 

x€[p,X]  xt[p,X] 

While  it  may  not  always  be  true  that  P (F  ,F. )  *  G(X-u)  for  the  situa- 

M  A 

tion  under  consideration,  we  are  able  to  state  the  following  result. 
Theorem.  Suppose  G  is  a  continuous  function  which  Is 
strictly  increasing  for  positive  values  of  its  argument. 

If,  for  each  positive  real  number  6  lens  than  a(q)-X  , 

G(X-n+6)  <  G(X-p)  +G(&)  ,  then  P (F  ,F.)  -  G(X-^>  . 

U  A. 

Proof.  If  x  e  (X,a(q)l  »  then  (x-X)  e  (0,a(q)-X]  ,  go  that 

x-X  may  be  taken  as  6  .  If  this  substitution  is  made  in  the  in¬ 
equality  on  G  ,  then  we  see  that  G(x-n)  <  G(X-u)  +  G(x-X)  ,  which 
is  equivalent  to  d(x;p,X)  <  G(X-^)  for  x  €  (X,a(q)]  ,  thus  com¬ 

pleting  the  proof  of  the  theorm. 

To  see  that  the  exponential  example  treated  in  equations  (1)  and 
(2)  is  indeed  covered  by  the  ineq’ality  of  the  theorem,  it  is  sufficient 
to  note  that  the  distribution  function  of  the  exponential  distribution 
given  by  (2)  satisfies  the  functional  equation 

G(x+y)  -  G(y)  -  G(x)[l-G(y)] 
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for  *11  x  >  0  and  y  >  0  .  Since  X  >  n  and  6  >  0  ,  we  tee  that  for 

the  G  of  equation  (2)  we  have 

G(X-n+6)  -  G(6)  -  G(>-m)[1-G(6)]  , 

and  hence  that  the  Inequality  of  the  theorem  Is  satisfied,  since  for 

any  6  >  0  ,  we  have  G(&)  >  0  • 
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IV.  TWO  THEOREMS  ON  SOLUTIONS  OP  DIFFERENTIAL-DIFFERENCE 


EQUATIONS  AND  APPLICATIONS  TO  EPIDEMIC  THEORY 


1.  Introduction 

We  present  two  theorems  that  provide  simple  Iterative  solutions 
of  special  systems  of  differential-difference  equations.  We  show  as 
examples  of  the  theorems  the  simple  stochastic  epidemic  (cf.  Bailey, 

1957,  p.  39,  and  Bailey,  1963)  and  the  general  stochastic  epidemic  (cf. 
Bailey,  1957;  Ganl,  1965;  and  Siskind,  1965),  in  each  of  which  we  let 
the  Initial  distribution  of  the  number  of  uninfected  susceptibles  and 
the  numb  ;r  of  infectlves  be  arbitrary  but  assume  the  total  population 
size  bounded.  In  all  of  the  references  cited  above  the  methods  of 
solution  Involve  solving  a  corresponding  partial  differential  equation, 
whereas  we  deal  directly  with  the  original  system  of  ordinary  dlf ferential - 
difference  equations.  Furthermore  In  the  cited  references  the  authors 
begin  at  time  t  =  0  with  a  population  having  a  fixed  number  of  unin¬ 
fected  susceptibles  and  a  fixed  number  of  Infectlves.  For  the  simple 
stochastic  epidemic  with  arbitrary  Initial  distribution  we  provide 
solutions  not  obtainable  by  the  results  given  by  Bailey  (1957  or  1963). 

For  the  general  stochastic  epidemic.  If  we  use  the  results  of  Ganl  or 
Siskind,  then  the  solution  of  the  problem  having  an  arbitrary  Initial 
distribution  would  involve  additional  steps  that  would  sum  proportionally- 
weighted  conditional  results. 
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Let  j(t)  end  x'(t)  denote  n  by  1  column  matrices  whose  i’th 

row  elements  are  the  real-valued  differentiable  function  x^(t)  and  its 

derivative  x^(t)  ,  respectively,  each  defined  for  all  t  >  0  .  Let 

initial  conditions  at  time  t  -  0  be  x(0)  m  a  ,  where  the  column 

matrix  a  has  as  i'th  row  element  the  real  number  a.  .  Let  B  be 

t*  i  r* 

an  n  by  n  triangular  matrix  whose  (i,J)'th  element  b(i,j)  ir  s 
constant  for  each  pair  (i,j)  and  in  particular  b(i,j)  *  0  for  i  <  j 

We  let  I  ■  (1 . n)  and  for  convenience  of  notation  we  denote  b(i,i) 

by  b,  for  i  €  I  .  Let  C  denote  an  n  by  n  matrix  with  c(i,j) 

as  (i,j)'th  element.  Occasionally  it  will  be  convenient  to  write 
b(i,j)  or  c(i,J)  as  bfi,J]  or  cf i , Jl  ,  respectively.  The  n  by 

1  column  matrix  with  i'th  row  element  expfb^t)  is  denoted  by  e(t) 

We  define  the  symbol  Bjtx)  to  be  equal  to  t  when  x  *  0  and  x"1 
when  x  ^  0  .  We  define  the  real-valued  function  52(x)  to  be  equal 

to  1  when  x  ■  0  and  x  *  when  x  4  0  •  Finally,  we  shall  make 

frequent  use  of  the  function  E(x)  defined  as  1  for  x  >  0  and  0 
for  x  <  0 


2.  Solutions  of  Some  Systems  of  Differential-Difference  Equations 

THEOREM  1.  Let  jf(0)  **  j  ,  and  for  t  >  0  let  x'(t)  =  ]}x(t)  ,  where 

(i)  b(i,j)  =  0  for  i-j  >  2  ,  and  (ii)  for  each  i  e  I  ,  b  *  b 

n  i  j 

for  at  most  one  J  e  In  and  J  4  1  •  Then  *(0  m  £e(t)  ,  where 


c(i.J) 


0  ,  i  <  j 

•i  •  i  -  J  -  1 

b(i,i-l)tc1(i-l,J)6l(bJ-b1)  -  c2(i-l,j)6^(bj-bi) 

+  c2(i-l,j)B1(bj-b1)t] 

•4  -  ^  cjd.u)  ,  i  -  j  >  i  , 

11*1 


(1) 

i  >  J 
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in  which  for  fixed  j  the  functions  c  and  c2  are  defined  recursively 
(in  i  )  as  the  term  independent  of  t  and  the  coefficient  of  t  , 
respectively,  in  c(i,j)  ;  i.e.. 


c(i,j)  -  c^i.j)  +  c2(i,j)t 


(2) 


(Thus,  in  particular,  c^(l,l)  =*  a^  ,  c2(l,l)  •»  0  ;  aud  for  i  >  1  , 

i-1 

Cjd.i)  *  ‘  c1(i,u)  ,  c2(i,i>  -  0  .) 

u=l 


blt 

Proof.  Note  that  x^(t)  *=  a^e  so  that  c(l,j)  is  equal  to  a^ 
when  j  =  1  and  0  when  j  >  1  .  For  fixed  integer  i^  where 

1  <  iQ  -  1  <  n  assume  equation  (1)  holds  for  positive  integers 
1  <  iQ  -  1  .  Consider  the  equation 


x  (t)  -  b  x  (t) 

0  0  l0 


b<Vlo 


V1 

1)  c(i0  -  1, j)e 

J-l 


V 


(3) 


The  term  in  which  j  •  Jq  5  i0  ”  *  00  t^ie  ri8bt  hand  side  of  equation  (1) 

bJ  1 
J0 

contributes  to  the  solution  of  J^(t)  *be  term  *(iQ«ig)e  •  We 

shall  show  that  Kdg.jg)  ■  c(iQ,j0)  .  Note  that 

bj  t  bj,  t  /-(b j  -  )t 

*<VVe  “e  °  b^10,10~  dt  ^ 

r  <bj  •  bi  )b  q 

+  c,<to-I.l0/«  0  0  It]  . 

Our  proof  distinguishes  three  cases:  Case  1.  b  ^  b  for 

J0  1 
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i  -  jQ  +  .  For  thl«  case  c2(iQ-l, jQ)  -  0  .  Therefore 

"<vv  -h(vio-i>'i<io-i'V<bj0  -  v'1  •  C*"J-  bi0-\  • 

Here  again  ^(lg-l,.^)  ”  °  *  but  now  *  b^0’ V^l^o’1'^*  ' 

Case  3.  b.  «  b  where  J.  <  k  <  K  •  Then 
J0  *0 

H‘,,i)-b«0.vi)it,<vi'V(b)o  •  v’1  ‘  c2<Io'I'Jo><\  ■  V 2 

} 

+  c2<I0‘1,J0)(bJo  -  bl0>'1‘1  • 

Mote  that  all  three  caaea  are  accounted  for  by  the  8^  symbol  as  used  in 
equation  (1)  where  t.,  >  jg  •  Therefore 


0  b  t  bi  * 

t  (O  -  V  c(i0,j)e  J  0  , 


"i  '  '  j_  '  0'J'  l 

°  J-l 

and  so  by  applying  the  initial  condition  x  (0)  ■*  a.  we  obtain 

*0  0 


*  _  »  -  '  c  (i  ,u)  .  This  completes  the  proof  of  Theorem  1. 

1  iQ  1'  O’ 

u»l 

THEOREM  2.  If  J(0)=a  ,  x'(t)-J}x(t)  for  t>0  ,  and  if  for  every 

pair  of  integers  a  <  0  such  that  bQ  -  bg  we  have  either  (i)  b(a  +  r  ><*) 
for  y  «*  1, . . .  ,3  -  a  ,  or  (ii)  b(g,  j)  »  0  for  j  ■  a . B*1  »  then 


$(t)  "  ££<0  where 


0  ,  i  <  J 

.  i  ■  J  ■  1 


c(i,J) 


j  I  X 

82(bj  *  ht)  y b(i,u)c(u, J)  ,  i  >  j 


a1  -  ^  c(i,u)  ,  i  *  J  >  1 
u-1 
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Proof.  Note  that  x^(t)  =  a^e  so  that  c(l,J)  is  equal  to  a^ 
when  j  <=  1  and  0  when  j  >  1  .  For  fixed  integer  i^  where 

1  <  1q  -  1  <  n  ,  assume  equation  (5)  holds  for  all  positive  integers 
1  <  iQ  -  1  and  let  R£  -  (J  :  bj  -  ,  1  <  j  <  iQ  -  1)  and 

R*  =  {J  ;  b  4  b  ,  1  <  J  <  ln  -  1)  •  Then 

0  J  0 


xj  (O  -  b  x  (t)  =  b(i  ,u)x  (t) 
in  in  1r>  Z-  0  u 


i  -1 

0  u  w  r 

=  b<in>u)  c(u,v)e  v 

L-j  ^  L  j 


u-1 

V=1 

(  r 

+ 

V1 

L  . 

VUrR* 

L  . 
vcR 

L  «s 

U=V 

b  t 

b(i  ,u)c(u,v)e  V  , 


Therefore 

r->  ^  b(in.  v)c(u,v)  b  t 

x  (t)  =  — TT — -  r — r-  e  V  +  D 

i0  Z  L..  <bv  hi  ) 

ME  R*  u=v  *0 

0 


/  r1  V"  b(Vu)c(u-v\  *-0 

+  '  -io  ■  L.  L.  <bv  -  bi  )  ^ 

u  v*  Rf  u«v  v  l0 


where 


V1 

r  r 

i  U 

L  > 

‘wR. 

u«v 

b(i„.u)c(u,v) 


b  t 

,v)jte  0 


If  R  is  empty,  then  the  term  D  does  not  appear  in  (6).  If 
l0 

R  is  not  empty,  then  for  fixed  j  c  R.  we  have 
10  0  *0 


L  b(i0>u)c(u, JQ)  -  0 


tf  either  b(iQtu)  -  0  for  u  *  J0,...,i0-1  ,  or  b(jQ  4  r.  jg)  =  0  for 

X  *  .  (The  latter  condition  is  seen  to  be  sufficient  when  we  use 

the  identity  for  i  >  J  : 

jB(bJ-bi)b(i,j)c(j,j)  ,  i=j4l 

c<i,J)  -  S 


6(bj-bi)fb(i,j)c(j,j)  4 


r 

u 

u=j4l 


b(i,u)c(u, j)l 


i  ■  ]  +  2 . n  .)  (8) 


Therefore  D  is  identically  equal  to  zero. 
It  is  now  easy  to  see  that  we  may  write 

\  b  t 

x.  (t)  =  c(i  ,v)e  v  , 

0  L\ 
v=l 

where 


0  ,  v  >  1, 


V1 


c(i.,v)  6-(b  -b  )  b(i  ,u)c(u,v) 

0  *  Y  1q  L  -  u 


u«v 


V1 

iv  -c(1°‘u>  '  v“0 

^  U  U— 1 


V  <  i 


and  so  the  theorem  is  proved. 


3.  Stochastic  Epidemics 

By  an  epidemic  population  we  shall  mean  a  well-defined  set  n  of 
elements  (individuals)  oj  defined  to  be  in  0  if  and  only  if  for  some 
time  t  >  0,fu  is  an  uninfected  susceptible  or  an  infective.  For  each 
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fi'  e  p  and  for  each  t  >  0  we  define: 


Wjfrst) 


W2(u>,t) 


w3(05,t) 


,  if  rjj  is  an  uninfected  susceptible  at  time  t 
,  otherwise 

.  if  fu  is  an  infective  at  time  t 
,  otherwise 

,  if  m  is  neither  an  uninfected  susceptible  nor 
an  infective  at  time  t 

,  otherwise. 


(9) 


We  shall  assume  that  the  number  of  elements  in  P  is  M  ,  a  finite 
positive  integer. 

Let 


r1 

R(t)  »  W  (<u,t) 

L  1 


S(t)  =  *Yo,t)  ^ 

L.  * 

coeO 


L(t)  a  W 

L  3 
fjoeP 


(10) 


Then 


M  -  R(t)  +  S(t)  +  L(t) 


(ID 
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We  shall  denote  the  size  of  the  epidemic  population  at  time  t  >  0  by 
N(t)  ,  which  consists  of  the  total  of  uninfected  susceptibles  and 
infectives  at  time  t  >  0  ;  i.e.. 

N(t)  -  R(t)  +  S(t)  .  (12) 

The  problem  we  consider  is  to  find  p  (t)  ,  the  probability  that 

ITS 

R(t)  -  r  and  S(t)  «=  8  ,  when  we  are  given  the  initial  distribution 

(pr#(0)}  and  information  about  the  infinitesimal  transition  probabilities 
for  an  to  move  amongst  the  three  states  of  being  (1)  an  uninfected 
susceptible,  or  (2)  an  infective,  or  (3)  neither  an  uninfected 
susceptible  nor  an  infective. 


4.  The  Simple  Stochastic  Epidemic 

In  the  simple  stochastic  epidemic,  which  has  been  extensively 
Investigated  by  Bailey  (1957,  1963),  there  is  a  positive  integer  N  such 
that  for  each  t  >  0  the  probability  is  one  that  N(t)  »  N(0)  N 
Therefore  S(t)  «•  N  -  R(t) 

When  we  make  the  usual  assumptions  (cf.  Bailey  (1957),  p.  39) 

★ 

about  the  infinitesimal  transition  probabilities,  then  we  obtain 

”;,H-r<*>  "  <r+l><K-r-DPr+l, *-.-!<*)  '  r<N'r>Pr ,N-r(t>  <13> 

for  r  m  0, 1, . . .  ,N  ,  where  pr> (t)  a  0  if  r  <  0  or  r  >  N  .  We 
write  the  initial  conditions  for  this  system  as 


Pr,H-r(0)  "  *r ,N-r  , 


r  m  0, 1 , . . . ,N  ,  where  each  a  >  0  and  a  ■  1  .  Thus 

r-0 


^Because  there  is  no  loss  of  generality  for  our  purpose,  we  have  assumed 
throughout  sections  4  and  5  that  the  infection  rate  is  equal  to  one. 
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In  addition  to  including  an  arbitrary  initial  distribution,  we  have 
introduced  for  completeness  the  case  in  which  r«N  ,  which  corresponds 
to  an  initial  population  completely  free  of  infectives. 

We  shall  now  put  this  problem  into  the  framework  of  the  theory 
of  section  2. 

LE>tiA  1.  For  each  ordered  pair  of  Integers  (r,N-r)  let 


where  k  £  k(r;N)  »  N-r+1  (16) 

Then  the  system  of  equations  (13)  with  initial  conditions  (14)  , 


where  p^  N_f(t)  £0  if  r<0  or  r  >  N  ,is  equivalent  to  the  system 


x^(t)  -  (N-k+2)(k-2)£(k-2)xk.l(i)  -  (N-k+1) (k*l)x^(t)  ,  (17) 

•  with  initial  conditions 


^(0)  -  ^  ,  (18) 

k  e  I  ,  where  n  =*  N  +  1 
n 

Proof.  If  we  make  the  indicated  change  of  variables  then  for  r  «*  0,1 . N 

Pr,N-r(*>  "  \<*> 

P;,H-r^>  “ 

Pr+l,N-r-l(t)  "  E<k-2)%-1<*)  ’ 

where  k  »  N  -  r  +  1  Thus  k  takes  values  l,2,...,N  +  lBn  ,  and 
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so  equations  (17)  and  (18)  follow. 


The  solution  of  the  simple  epidemic  is  given  by 
THEOREM  3.  If  we  ha^e  the  system  of  equations  (17)  with  initial  conditions 
(18),  V.  €  I  ,  then  x(t)  =  Ce(t)  ,  where  b^,  =  -(N-i+1)  (i-1)  and 
c(i,J)  is  given  by  eauation  (1),  in  which  b(i,i-l)  =  (N-i+2) (i-2)E(i-2) 
Proof.  We  need  only  note  that  for  i  =  l,...,n  ,  b^  =  bN_3+2  and 

b^  A  bj  for  j  *  N-i+2  ,  j  4  i  .  Thus  Theorem  1  applies. 

Example .  We  illustrate  the  simplicity  of  the  theory  by  showing  the 
details  of  the  example  in  which  N  =  6  ,  and  initial  distribution 

(agg»  *25’  *24 *  *33’  *42’  *51’  *60^  =  *30,  *25,  .15,  .10,  .10)  . 

If  we  use  Lemma  1,  then 


xj(+) 


x’(t) 


x^(t) 


*6(*> 

x^(t) 


j 

If 

x1(t)' 

al' 

i 

i 

0-5  !  | 

i  I 

X2(t)j 

a2 

\ 

j 

5  -8  ; 

x3(t) 

a3 

1 

i 

f 

;  „ 

i 

8  -9  i 

j 

i 

x^(t) !  and 

a4 

i 

| 

i 

9  -8 

| 

x5(t)| 

°5 

i 

8  -5  1 

x6(t) 

a6 

. 

5  0j; 

x?  (-t)j 

“7. 

By  applying  Theorem  3  we  obtain: 


c(l,l)  *  Jq  and  c(i,l)  =  0  for  i  >  1 


c(2,2)  *»  1Q 

c(3,2)  -5(^)4) 

c<4,2)  -  8(|)<|)  =  \ 

c(5,2)  -  9(3X3)  =  1 

c(6,2)  *  8(1)61(0)  *  8t 

2 

c (7 , 2)  =  -  5(8) (-  3)  +  5(8) (-  3>c  =  -8t 

=».3)  » is  -  ;  -  -  is 

c<4,3)  .  »(-  -|o> <T>  -  -  15 
c(5,3)  -  9(*|3)61(0)  =  '3  t 

c(6,3)  =  -  8<*|)(‘±)  +  8(‘|)(*^)t  =  yf  +  ^  c 

c (7 , 3)  =  5(|§)(-  3)  -  5(^f)(*|)  +  5(^|)(-  |)t  =  -  2t 

c(4  4)  -  i  -  i  +  —  =  — 

CV '  ’  4  3  15  20 

c(5.4)  ~  9(|q)(-  1)  =  -  33 
c(6,4)  =  8(-  33)  (-  \)  =fo 
c(7,4)  =  5(|q)(-  3)  =  -  3 

c(5,5)  =  Jo  '  1  +  20  "  '  4 
c(6,5)  =  8(-  £)(-  ±)  =3 
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c<7,5) 

c<6,6) 

c<7,6) 

c(7,7) 


5(|)(-  J> 


5_ 

12 


10 


16 

15 


9_ 

10 


5(-  h 

30m  5' 


2  _76 

3  "  ‘  30 

76 

30 


8  11  1  5_  76  9_ 
5  12  2  +  12  ‘  30  “  10 


Therefore  Che  solution  is 


p60 

«, 

!  6 

1 

P51 

X2 

0 

6 

1 

P42 

X3 

0 

10 

-1 

P3, 

*4 

1 

"  60 

0 

20 

-8 

3 

P24 

*5 

0 

60 

-72t 

-27 

-15 

P15 

X6 

0 

480t 

64+192t 

54 

40 

-152 

P06. 

V 

.  0 

-96-480C 

-55-120t 

-30 

-25 

152  54  / 

6 

1 

0  6 

e^ 

0  10  -1 

e-8* 

0  20  -83 

-9* 

e 

0  60  -15  -27  -72 

te*8* 

0  -152  104  54  192  480 

te'5* 

.0  56  -80  -30  120  -480  54. 

.  1  - 

Mote  thst  of  the  28  entries  calculsted  for  CJ  only  5  hsd  4  0 


1 


■5t 


-8t 


-91 


•8t 


-5t 


1 
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Once  Che  individual  probabilities  are  known  Chen  it  ia  an  easy 
matter  Co  obtain  ocher  derived  quantities  such  as  the  duration  of  the 
epidemic,  the  mean  number  of  infectives  at  time  t  ,  and  the  distri¬ 
bution  of  the  total  size  of  the  epidemic.  For  example,  the  distribution 
function  of  time  to  extinction  may  be  read  directly  from  the  above 
solution  as  G(T)  «  P0^(T)  B  *7^)  *or  T  >  0  »  •*»«*  G(T)  -  0  for 


T  <  0  . 

5.  The  General  Stochastic  fpldemle 

In  the  general  stochastic  epidemic,  which  has  recently  been 

investigated  by  Gani  (1965)  and  Siskind  (1965),  the  probability  is  one 

that  for  tj>tj>0  .  Nftj)  <  N(t^)  .  We  shall  assume  the  initial 

population  bounded  in  the  sense  that  there  exists  a  positive  integer  N 

such  that  the  probability  is  one  that  N(0)  <  H  .  Let  S'  be  the  set 

*  II 

of  ordered  pairs  of  Integers  {(r,s)  :  r  >  0  ,  a  >  0  ,  r  +  s  <  M) 

If  we  make  the  usual  assumptions  (cf.  Bailey  (1957),  p.  53)  about  the 
infinitesimal  transition  probabilities  then  we  obtain 

pra(t)  *  <r+1)<*'1)prfl,s-l(t*  "  +  P^DPr  ^d)  (20) 


where  p  is  the  removal  rate  and  p  (t)  e  0  if 

rra' 

write  the  initial  conditions  for  this  system  as 


(r,«)  $ 


We 


P  (0)  -  a 

rs  rs 


(r,s)  e  ,  where  a^  >  0  and  ^ 


a  -  1 
rs 


(r,a)eS^ 


(21) 


In  order  to  put  this  epidemic  problem  into  the  framework  of  the  theory 
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of  section  2  we  construct  s  counting  mechenlsm  for  the  equations  and 
the  variables.  This  Is  done  in  the  following  lemma  whose  proof  Is  left 
to  the  reader. 

LEMMA  2.  For  each  non-negative  Integer  N  ,  let  denote  the 

set  of  ordered  triplets  of  Integers  (k,r,s)  ,  where  k  >  0  ,  r  >  0  , 
s>0,r  +  s<N  ,  and 

k  £  k(r , s ;N)  -  (N+l)(N+2)/2  -  (N-fl)r-  s  +  (r-l)r/2  .  (22) 


Then  contains  exactly  n  *  (N+l)(N+2)/2  ordered  triplets  and  for 

each  positive  integer  k  <  n  ,  there  exi«r«  one  and  only  one  ordered 
pair  of  non-negative  integers  (r,s)  such  that  (k,r,s)  e  Su  . 

Therefore  for  each  pair  (r,s)  ?  S'  one  can  find  k(r,s;N) 

N 

It  sight  be  worthwhile  to  point  out  that  the  converse  problem  can 
also  be  neatly  treated;  namely,  for  each  positive  integer  k  <  n  , 
let  u  be  the  greatest  Integer  which  is  less  than  (1  +  V8k+l)/2 
Then  r  =»  N  +  1  -  u  and  s  ■  u(u  +  l)/2  -  k  .  (Later  in  the  state¬ 
ment  and  proof  of  Theorem  4  it  will  be  convenient  to  use  the  notation 
(r^.a^)  in  order  to  indicate  the  one-to-one  correspondence  between 
k  £  k(r,s;N)  and  (tk»8k)  •) 

We  are  now  prepared  to  effect  a  change  of  notation,  which  we  do 
in  the  following  lemma. 

LEMKA  3.  For  each  ordered  pair  of  integers  (r,s)  let 


Prg(t)  ■  £(r)£(s)£(N-r-s)xk(t) 


(23) 
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where  k  s  k(r,a;N)  is  given  by  (22).  Then  the  ayatem  of  equation*  (20) 
with  initial  conditioua  (21),  where  Pra(fc)  *0  if  (r,a)  i  ,  ia 
equivalent  to  the  system 

^(t)  +  ■(rH>)xk(e)  -  (r+l)(a-l)£(a-l)xk_N+r(t)  +  P(a-H)£(N-r-l-a)xk_1  (t)  (24) 

with  initial  conditions 


^(0)  -  «k  .  (25) 

k  €  I  ,  where  n  -  (N  +  1 ) (N  +  2)/2  and  (k,r,a)  €  S„  . 
n  n 

Proof.  If  we  make  the  indicated  change  of  variablea  then  for  each 
(k,r, a)  c  S„  we  obtain 


Pr.(t>  -  v*) 

Pr.(t)  -<(*) 

pr+l,a-l(t)  “  £<*‘l>Vl»r(e) 
pr,.+l(t)  *  €(N-r-l-a)Vl(t) 


(26) 


and  so  the  concluaion  of  the  lemna  follows  turned  lately. 

We  now  state  the  solution  of  the  general  stochastic  epidemic  i~ 

THEOREM  4.  If  we  have  the  system  of  equations  (24)  with  Initial 

conditions  (25),  k  €  I  ,  where  n  ■  (N  +  1) (N  +  2)/2  and  (k,r,s)  c  Su  , 

n  N 

and  if  P  is  such  that  for  (r,s)  e  ,  (r’,s')  €  ,  s  i  0 

and  s'  0  ,  we  have  a(r+P)  -  *'(r'4P)  only  if  a  ■  s'  ,  then 

x(t)  -  Ce(t)  ,  where  for  any  T  e  I  ,  b  -  -  a„(r  +  P)  ,  and 

n  i  f  [ 
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0  .  k  <  J 


*1  • 


k  -  J  -  1 


e(k, j)  -  ■{  [ (r+l)(a-l)E(e-l)c(k-lf+t,))£(k-IWr-j) 

+  p(a+l)8(N-r-l-*)c(k-l, J)]62(bj  -  bfc)  ,  k  >J 


(27) 


a  -  C(j,u)  ,  k  -  j  >  1  . 

u-1 


Proof.  If  (r ,a)«  and  a  >  0  ,  then  bk^r>s;Hj  ■  "  •(*  +  p>  *re 

all  dlatlnct.  If  a  -  0  ,  then  for  r  -  0,1 . N  we  have  bk(rt0;Nj  *  0  • 

The  following  argument  holda  for  each  poaitlve  integer  r  <  N  .  In 
equation  (24)  the  only  poaaible  noo-xero  coefficlenta  of  ^^q.h)  ar® 

b[k',k(r,0;H)]  and  b[k",k(r,0;N)l  ,  where  k'  -  k(r,0;N)  +  1  and 

k"  -  k(r,0;N)  +  H  -  rk„  ,  We  aee  that  k'  -  k(r-l,H-r+l;N)  and 
k"  -  k(r-l,l;H)  .  Therefore  b[k',k(r,0;N)]  «  Pt (N-r+l)+l]E(N-(r-l)-l-(N-r+l)) 
P(H-r+2)£(-l)  ■  0  and  b(k",k (r,0;N)]  -  [ (r-l)+l] (l-l)E(l-l)  -  r(0)£(0)  -  0  . 

Thua  condition  (i)  of  Theorem  2  la  aatiafied.  Pinally,  by  applying 
equation  (5)  we  get  the  conclualon  of  our  preaent  theorem.  Note  the 
reaemblanea  between  equation  (24)  and  equation  (27)  for  the  caae  k  >  J  . 

The  queation  arlaea  aa  to  whether  or  not  there  la  anything  dlatinctlve 
about  the  choice  of  If  ,  the  bound  on  the  initial  total  population  aixe. 

The  following  ea^odding  theorem  anawera  thia  queation  in  the  negative. 

THBOREM  5*  the  gw  rsal  etc'chaatic  epidemic  with  total  initial  population 
else  bounded  by  K.  <  cf„  may  be  treated  aa  one  with  initial  total  pop- 
u let ion  aixe  boused  by  Kj  and  initial  condltiona  aatiafying 
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a  f  x ”  0  for  all  non-negative  Integer  pair*  (r,a)  such  ♦•hat 

ic\r  |  •  * 

r  +  a  >  Hl  . 

Ve  leave  the  details  of  the  proof  to  the  Interested  resder. 

Essentially  what  must  be  shown  Is  that  for  each  non-negative  Integer 

pair  (r,a)  such  that  r  +  s  <  N.  we  have  x.  .  At)  *» 

“  *  *ir,s;H^ 

x.  /  „  .  (t)  for  all  t  >  0  .  Also  for  each  non-negative  Integer 

*ir*s»«2'  ” 

pair  (r,s)  such  that  N.  <  r  +  a  <  N,  we  have  .(t)  £  0 

l  —  l 

The  special  cases  considered  by  Canl  and  Siskind  are  included 
in  Theorem  A  by  simply  choosing  the  initial  conditions  appropriately. 

In  fact  if  a  computing  program  Is  worked  out  for  fixed  N  ,  then  any  of 
their  casws  in  which  R(0)  =  r^  and  S(0)  »  sQ  ,  where  r^  +  sQ  »  N  , 


may  be  obtained  by  letting  p  (0)  »  a 

Vo  ro8o 


and  setting  all  other 


a^'s  equal  to  zero.  Furthermore,  by  Theorem  5,  any  case  in  which  the 

probability  is  zero  that  the  sum  of  R(0)  and  S(0)  exceeds 

may  be  obtained  from  the  program  written  for  Nj  by  suppressing  from 


the  program  all  terms  involving  triplets  (k,r,s)  g  SH 

Vi 

r  +  s  >  N, 


such  that 


’1  * 

Conversely,  if  one  wishes  to  use  the  results  of  Canl  or  Siskind 
to  solve  a  general  stochastic  epidemic  with  arbitrary  initial  distri¬ 
bution  then  one  can  obtain  prg(t)  by  evaluating 


(r0’V6SN 


PfR(t)  -  r,S(t)  -  slR(O)  -  rn,S(0)  -  >Ja 

U  0  Vo 
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where  P[A|B]  is  the  conditional  probability  of  A  given  B  .  This, 
of  course,  would  be  considerably  store  difficult  than  if  one  were  to  use 
our  Theorem  4  directly. 

Example.  We  illustrate  the  simplicity  of  the  theory  by  showing  the 
details  of  the  example  in  which  N  m  2  ,  p  -  2  ,  and  initial  distri¬ 
bution  (*qq»  *q2'  *10’  *11*  *20^  m  *20,  .30,  .25,  .15,  .10)  . 

Then  SN  -  {(1,2,0),  (2,1,1),  (3,1,0),  (4,0,2),  (5,0,1),  (6,0,0))  . 

If  we  use  Ltwms  3  then 


By  applying  equation  (27)  we  obtain 


c(l.l) 

c(2,l)  -  (  0  +  0c(l,l))/3  -  0 
c(3,l)  -  (  0  +  2c(2,i>62(0)  -  0 
c(4,l)  -  (lc(2,l)  +  0  )/4  -0 
c(5,l)  -  (0c(3,l)  +  4c(4,l J/2  »  0 
c(6,l)  -  (  0  +  2c(5,l|R2(0)  -  0 
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c(2,2)  ”  20  ”  c^2,1^  “  20 

c(3,2)  -  (  0  +  2c(2,2))/(-3)  -"jq 

c(4,2)  m  (lc(2,2)  +  0  )/l  “  20 

c(5,2)  -  (0c(3,2)  +  4c(4,2))/(-l) 
c(6,2)  -  (  0  +  2c<5,2))/<-3)  -  | 

c(3,3)  ■  -  c(3,l)  “  c(3,2)  “  20 

c (4 , 3)  -  (  0  +0  )/(4)  -  0 

c(5,3)  -  (0c(3,3)  +  4c(4,3))/<2)  -  0 

c(6,3)  «  (  0  +  2c(5,3))&2(0>  =  0 

c(4,4)  ■  Yq  -  c(4,l)  -  c(4,2)  -  c<4,3)  *20 
c(5,4)  -  (  0  +  4c(4,4))/(-2)  -  |o 

c(6,4)  «  (  0  +  2c(5,4))/(-4)  -  |o 

c(5,5)  -  \  -  c(5,l)  -  c(5,2)  -  c(5,3)  -  c(5,4)  »  ^0 
c(6,5)  -  (  0  +  2c<5,5))/(-2)  -  - 

5 

c(6,6)  -  -  2  c(6,u)  =  jo 

U=1 
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Therefore  Che  solution  is 


p20 

*1 

> 

2 

r  1  1 

1  -3t 

P11 

X2 

0  3 

e 

p10 

s 

X3 

1 

ol 

0-2  7 

1 

-4t 

P02 

x4 

0  3  0  3 

e 

-2t 

P01 

X5 

0  -12  0  -6  22 

e 

P00/ 

X5 

0  8  0  3  -22  11  J 

-  1  - 

The  same  comment*  which  were  made  in  the  example  of  section 
4  regarding  the  ease  of  finding  other  derived  quantities  from  the  above 
solution  apply  here.  For  example,  the  distribution  function  of  time 
to  extinction  obtained  from  the  above  solution  is  G(T)  -  0  for  T  <  0 
and  for  T  >  0  , 

2 

C(T)  -  Y  PrO(T)  "  1  +  +  3e"4T  '  22e'2T>  * 

r-0 
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